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Abstract 

Nonlinear Schrodinger / Gross-Pitaevskii equations play a central role in the understanding of nonlinear 
optical and macroscopic quantum systems. The large time dynamics of such systems is governed by 
interactions of the nonlinear ground state manifold, discrete neutral modes ("excited states") and dispersive 
radiation. Systems with symmetry, in spatial dimensions larger than one, typically have degenerate neutral 
modes. Thus, we study the large time dynamics of systems with degenerate neutral modes. This requires 
a new normal form ( nonlinear matrix Fermi Golden Rule) governing the system's large time asymptotic 
relaxation to the ground state (soliton) manifold. 
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1 Introduction 



Nonlinear Schrodinger / Gross-Pitaevskii (NLS/GP) equations are a class of dispersive Hamiltonian partial 
differential equations (PDEs) of the form: 

idMx,t) = -A^(x,t) + {V{x) - f{\i^{x,t)\^))i^{x,t). (1.1) 

Here, ip = tp{x, t) is a scalar complex-valued function of position, x G M*^ and time, t G M. The function 
y : R'^ ^ R denotes a linear potential and / : R+ — > R, a nonlinear potential. For example, V can be taken 
to be a smooth, non-positive potential well, with rapid decay as \x\ — > oo and = —g g constant. 

For 5 > 0, the nonlinearity is called repulsive or defocusing. For 5 < it is called attractive or focusing. In 
this paper, we focus on spatial dimensions d > 3. Precise hypotheses on V and / are given below. We are 
interested in the initial value problem (IVP) for (11. ip with finite energy data V'(a^)O) and solution 'ip{x,t), 
which are sufficiently regular and decaying to zero as \x\ 00. A precise well-posedness result is cited 
below; see Theorem 13.11 

NLS /GP equations play a central role in the understanding of nonlinear optical [4H 21 [32] and macro- 
scopic quantum systems; see, for example, |19j . A striking and important feature of NLS/GP is that it 
can have localized standing waves or nonlinear bound state solutions, some of which are stable and play a 
central role in the general dynamics. In particular, for a wide variety of potentials and nonlinearities there 
exists an interval T C R such that for any A G T, (jl.ip has nonlinear ground state solutions. These are 
solutions of the form 

ip{x,t) = e'^*^^{x), 

where 

-A^^ + V^^-f{\cl)^f)cl)^ = -X^^ (1.2) 

with 4>-^ E and <j)-^ > 0. 

The gauge (phase-translation) invariance of ()1.1|) . 

tp e'-'^ip, 7 e [0, 2tt), 

generates a nonlinear ground state or "soliton" manifold: 

Mj := {e'^0\ A e T, 7 G [0, 2^)}. (1.3) 

If V is identically zero, then NLS /GP admits a larger group of symmetries and the definition of soliton 
manifold (which exists in the focusing case, g < 0) is naturally extended to incorporate these additional 
symmetries; see, for example, [571 [26], 

Orbital Stability: The soliton manifold Aij is said to be orbitally stable if any initial condition ipQ, which 
is close to Aij in H^, gives rise to a solution tp{t), which is close for t ^ 0. There is an extensive litera- 
ture on the orbital stability of the soliton manifold. For the case, V = 0, orbital stability (stability modulo 
spatial and phase translations) of global energy minimizers was proved in [9j by compactness arguments. 
In [56l I57j it is shown that positive solutions, which are index one critical points (Hessian with one strictly 
negative eigenvalue) and satisfying the slope condition H: 

/ \cp\x)\^ dx > 0, (1.4) 
"A J^d 

^The term soliton sometimes refers, more specifically, to particle-like solutions of completely integrable PDEs. 

= —A is the typical definition of soliton frequency. Therefore the slope condition, (|1.4|l . often appears as a rate of change 
with respect to /i being negative. 
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are orbitally stable. For V = and /(IV'P) = —g\ip\'^T9 < (focusing case), (|1.4p is equivalent to 
a < 2/d. Orbital stability of solitary waves of NLS/GP for a class of potentials, V, was studied in [39] and 
in a semi-classical setting by J35]. A general formulation of a stability / instability theory is presented in [26]. 

Asymptotic stability: We say the soliton manifold, Aij is asymptotically stable if ipo close to Aij in 
a suitable norm implies that 'ip{t) remains close to and converges to Mj (in a possibly different norm), as 
t tends to infinity. 

Are solitary waves asymptotically stable? This is a local variant of the problem of asymptotic resolution^U\, 
i.e. whether general initial conditions resolve into stable nonlinear bound states of the system plus dis- 
persive radiation. A great deal of progress has been made on this problem in recent years. The study 
of asymptotic stability of solitary waves was initiated in [H] ; see also [51 [Ml [2H1 [55] . In the translation 
invariant case, asymptotic stability was then investigated by [6]. Asymptotic stability analysis requires 
two new analytical features: one dynamical systems and the other harmonic analysis / spectral theoretic. 

First, since we do not know in advance which nonlinear ground state in Mx emerges in the large time limit, 
a decomposition with flexibility allowing for the asymptotic soliton to dynamically emerge is requirecH. 
To this end, the solution is decomposed in terms of a motion along the soliton manifold and components 
(symplectic or bi-) orthogonal to it. Dynamics along the soliton manifold, A4i, are governed by modulation 
equations] see, for example, [56 l [211 [31] . 

Secondly, in order to prove convergence to the soliton manifold, Mj, we need to show that the devia- 
tion of the solution from Mj decays with advancing time. This requires time-decay estimates [U'{W^) or 
space-time norms) for the linearized (about the soliton) propagator on the subspace (symplectic or bi-) 
orthogonal to the discrete spectral subspace. The discrete subspace is the union of a zero frequency mode 
subspace spanned by infinitesimal generators of the NLS/GP symmetries (translation, gauge) acting on (^•^, 
and often a subspace of neutral modes (sometimes called internal modes) with non-zero frequencies. 

Since a typical perturbation of the ground state solitary wave in Mx excites all discrete spectral 
components, one must understand the mechanisms, due to which these do not interfere with the asymptotic 
convergence of ip{x, t) to ^Ax■ In brief: Concerning the zero modes, the choice of modulation equations 
"quotients out" the zero modes; perturbations exciting these, induce motion along the soliton manifold. 
And concerning the non-zero frequency neutral modes, these are shown to damp to zero, as t ^ oo, due 
resonant nonlinear coupling of discrete to radiation modes. Related to this is a further dynamical systems 
aspect of the analysis. The neutral mode amplitudes are governed by nonlinear oscillator equations, coupled 
to a dispersive wave field. Near-identity changes of variables are used to put the system in an appropriate 
normal form, wherein the mechanism of energy transfer from the neutral modes to the evolving soliton 
and propagating radiation, is made explicit. Energy transfer shows up as an explicit (nonlinear) damping 
term in the normal form; see the discussion below. The positive damping coefficient (matrix, in the present 
work) is a nonlinear variant of Fermi's Golden Rule [llj . See [6], regarding the dynamics near solitary 
waves of the translation invariant NLS equations and [47] for "breathers" of a class of nonlinear wave 
equations. In [l8j this mechanism was proved to be responsible for ground state selection in NLS/GP 
equations; see also [58]. Experimental verification of the prediction in [48[ [49] is reported in [40]. Related 
work on resonant radiation damping appears in [541 [71 [531 [5T1 [161 [T7] . The role of the Fermi Golden rule in 
the non-persistence of coherent structures for nonlinear wave equations was first demonstrated, via Floquet 
analysis, in There is a close relation to the perturbation theory of embedded eigenvalues for linear 

problems [371 [Ml [l8] . 

■^The case of integrable systems, such as Id NLS {V = 0, /(iV'P)'0 = I'^/'P^) is an important class for which it is possible 
to determine the emerging coherent structures from the scattering transform of the initial data. 
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The above works on nonlinear resonance required that the neutral modes frequencies (a) lie sufficiently 
close to the essential spectrum and (b) be of geometric multiplicity one. For example: for the cubic non- 
linearity, /(IV'P) = "S'lV'P) close means: coupling of to radiation modes occurs at order l^p. The situation 
where simple neutral modes with a large spectral gap has been studied in [25l |23l |22l [iTl E] • Here, coupling 
of the discrete to continuum modes occurs at some high order in g. Thus, the normal form expansion gives 
a damping term at some even order: Ig]"^^, k > 2. 

Results of this paper - systems w^ith degenerate neutral modes: 

An important situation, not covered by previous results, is the dynamics in the presence of degenerate 
neutral modes. This case arises naturally in systems of spatial dimension d > 2 with symmetry. For 
example, if the potential is spherically symmetric, V = V{\x\), then the first and higher excited states are 
degenerate, the degree of degeneracy related to the order of the associated spherical harmonics. Another 
interesting class of examples is a class of double-well potentials; see Appendix 113. 1[ 

In this paper we prove the asymptotic stability of the ground state / soliton manifold, of NLS/GP 
when the linearized spectrum has degenerate neutral modes. We show that the solution has three inter- 
acting parts: 

(i) a modulating soliton, parametrized by motion along Mx 

(ii) oscillatory, spatially localized, neutral modes which decay with time and 

(iii) a dispersive part, which decays in a local energy norm. 

The neutral modes and dispersive waves decay via transfer their mass to the soliton manifold or to spa- 
tially infinity. Additionally, degenerate neutral modes are coupled and exchange mass among themselves 
in addition to the soliton and radiation. These degenerate modes cannot be viewed as very weakly coupled 
"oscillators" |51) . We require instead a new normal form expansion. This is related to ideas developed in 
[34], where a parametrically forced linear Hamiltonian PDF was considered, and a normal form, uniform 
in discrete eigenvalue spacing, was required. 

We outline the perspective we take and give a rough form of the main theorem, Theorem 17. 11 Consider 
NLS/GP, where —A -|- V has a ground state, Coi^) > 0, eo < and a degenerate excited state, whose 
energy, ei, cq < ei < 0, is assumed sufficiently near the zero. Typical solutions of the linear Schrodinger 
equation, evolving from localized initial data V'Oi = exp[— z(— A -|- V)t]TpQ, will evolve a (typically) 

time quasi-periodic superposition of spatially localized ground and excited time-periodic states, plus a part 
which disperses to zero, i.e. tends to zero as t advances in Lf^^. This picture emerges from the spectral 
decomposition of —A + VinL'^, with respect to which the bound state projections of the solution evolve 
as independent oscillators and the continuous spectral part of the solution has a character, qualitatively 
like a solution to the free Schrodinger equation. 

For NLS/GP, e.g. —g\ip\'^ip, g ^ 0, the dynamics of discrete and continuum modes are coupled. We 
consider an appropriate open set of initial conditions near the soliton manifold. In contrast to the 
linear Schrodinger equation we show that the solution converges to a nonlinear ground state. To see 
this, we view NLS/GP as a infinite dimensional Hamiltonian system comprised of two subsystems: (i) a 
finite dimensional system governing dynamics along the soliton manifold, Mj, parametrized by {X{t),^{t)), 
the zero modes amplitudes (01,02), and the neutral mode amplitudes, z = {zi,Z2,- ■ ■,Zn)'^, and (ii) an 
infinite dimensional dispersive Schrodinger wave equation. A very detailed analysis of this coupled system 
(the bulk of this paper) yields the following (rough) form for the asymptotic behavior of small amplitude 
solutions of NLS/GP: 

'Main Theorem' 1.1. Consider the initial value problem for NLS/GP. Assume arbitrary localized initial 
data, which are sufficiently near a small amplitude nonlinear bound state, . Then the solution of 
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NLS/GP evolves as a modulated soliton plus decaying error having the following form: 

^{t) = e*/o^W'^^e*W*)+''2(^W,2(i))] X [ ^Ht)+aiizit),zit)) + 0{\z{t)\) + R{t) ], (1.5) 

Here, X{t) Aqo- represents a localized, nonspreading but decaying part, satisfying \z{t)\ < 

C{t)~^^'^. Also, Oj = aj{z,z) = 0{\z\'^). R{t) represents a spreading, dispersively decaying part, and tends 
to zero as t ^ oo in Lf^^; more precisely, \\{x)~'^ R{t)\\2 — > 0, v > 0. 

For the precise statement, see Theorem 17.11 

A key part of the proof of Theorem lT.ll is to show that \z{t)\ tends to zero and that X{t) has a limiting 
value Aoo G T as f tends to infinity. We prove the latter, by showing dtX{t) S L-^(M"'"). We have two 
comments on the approach of this article to these issues. 

(1) New normal form: We show that there exist a non-negative symmetric matrix r{z,z) = 0(|2:p) 
and a skew symmetric matrix A{z,z) = 0(|zp) (see (17. 3p below) such that 

dtz = -iE{X)z-r{z,z)z + A{z,z)z + o(^{l + t)-^-^^ , S>0 (1.6) 

The matrix F is defined in terms of the spectral decomposition of the L{X) = JH(X), the generator of the 
linearized flow about the nonlinear bound state, (j)^; see section [5j Our analysis requires that F = F(z, z; X) 
is positive definite for an open interval of A— values. A variant of this hypothesis appears in previous work 
[Ml ISl El [53l EU EH EH [Ml [17] it is expected to hold, in some sense, generically. In section [6| we state 
hypotheses under which positive definiteness holds for class of potentials of double-well type, constructed 
in section llS.ll This hypothesis, denoted FGR ( see (16. 7116. 9p ), is a nonlinear variant of the Fermi Golden 
Rule [371 [IS] • We note that for finite dimensional Hamiltonian systems a damping term is absent ; it 
would violate phase-volume conservation. This term arises due to nonlinearity induced coupling between 
discrete and continuous spectral (radiation) modes, a phenomenon associated with continuous spectra, 
arising in PDEs on spatially infinite domains; see [471 158j . We show that ()1.6p and FGR (see (16. 7116. 9p 
) , imply the bound \z{t)\ = 0{t^2'j, For the case of multiple simple bound states with well-separated 
frequencies, a system of type (jl.6p holds with F, a diagonal matrix [51]. Equation (jl.6p can be viewed as 
a new normal form, a special case of one valid uniformly in neutral mode eigenfrequency-separation 

(2) Choice of basis for the neutral mode subspace: We prove that X{t) approaches some Aoo 

t oo, by proving that dtX{t) is integrable. If there are n simple well-separated neutral modes, one 
initially finds 

n 

dtX{t) = ^a„|z„|2 + 0(t-i). 

Since we expect \zm\ = 0{t^'i) we can not conclude integrability of dtX{t). However, it can be shown 
that after near identity change of variables: z ^-^ z + 0(|zp), we can take Om = 0; see the normal form 
expansion in [25l |23l [48] . In the degenerate (similarly not well-separated) X{t) satisfies the equation 

dtX{t) = ^am,kZmZk + 0{t~^). 

m,k 

In the present paper we show, very generally, by appropriate choice of neutral subspace basis we can take 
am,k = 0. 
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Finally, we expect that our techniques can be extended to more complicated situtations, e.g. where 
coupling of some neutral modes occurs at higher order in the nonlinearity. 

Outline of the paper: The paper is organized as follows. In Section [2] display notation, which is 
often used. Section [3] is a brief section outlining structural properties of NLS / GP and gives a statement 
of a basic well-posedness result. Section |1] introduces solitary waves (solitons) in the regime of weak nonlin- 
earity. Section [5] has a detailed discussion of the spectral properties of L{X) = JH{X), the generator of the 
linearized dynamics about the soliton: zero energy subspace, degenerate neutral subspaces and continuous 
spectral subspaces. Projections associated with theses subspaces are defined and decay estimates of the 
linearized evolution on the continuous spectral subspace are recalled. In section [6] the Fermi Golden Rule 
matrix, F, is introduced explicitly, (j6.8p . The detailed calculations, proving symmetry and non-negativity 
are given in the appendix; section [13.21 The main theorem, requires positive definiteness of F. Proposition 
16.21 is a result reducing the required positive definiteness to a condition involving the spectral properties 
of —A -|- V. Section [7] contains a statement of the main theorem, Theorem 17.11 In section [8] we give a 
more precise formulation of Theorem 17.11 This formulation makes explicit the dynamical (modulation) 
equations for the solitary wave parameters, the neutral mode amplitudes and the dispersive part. These 
are proved via normal form methods in sections [9] and [TOj In section [11] we prove the main Theorem 17.11 
in the setting of reformulated Theorem 18.11 Section [13] contains some important calculations used in the 
body of the paper. Of particular interest is the appendix of section 113. 1[ where a class of double-well 
three-dimensional potentials is constructed, to which we apply Theorem 17.11 
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2 Notation 

(1) a+ = max{a, 0}, [r] = maxfg^ {f < r} 

(2) 3^2 = real part of z, = imaginary part of z 

(3) Multi-indices 

w = {wi,...,wn) G C^, w = {wi,...,wn) 
a e N^, = • • • 4^ 
|a| = |ai| + ... + |ajv| 

z denotes the vector of neutral mode amplitudes 

^ denotes the vector, whose j*'* entry is the j^^ neutral vector-mode of JL{\), 

(4) Qm,n denotes an expression of the form 

TV 

Qm,n = Yl 1a,bZ''z^= ^ Qa,6 H ^fc' 
|a|=m, \b\=n \a\=m, \b\=n k=l 



(5) 



(6) aess{L) = (Tc{L) is the essential (continuous) spectrum of L, 
o'disciL) is the discrete spectrum of L. 

(7) Riesz projections: Pdisc{L) and Pc{L) = I - Pdisc{L) 
Pdisc{L) projects onto the discrete spectral part of L 
Pc{L) projects onto the continuous spectral part of L 

(8) 

{f,9) = J fix) 9{x) dx 

(9) 

11/11^ = / \f{xWdx, l<p<o 

(10) 

WfWns. = [ (x)^(/-A)f/(x)' 



dx 
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3 Hamiltonian Structure 

NLS/GP can be expressed as a Hamiltonian system: 

where the Hamiltonian energy, £\-]^ is defined by 

£[,P] = £[^M= j QvV-VV; + \v{x)^i, - F(V'V^)) dx, 

Equation ( II. ip is a Hamiltonian system on Sobolev space H^(M.'^, C) viewed as a real space H^{M.'^, M) © 
Fi(R'^,IR), i.e. H'^iM.'^X) 3 / ^ (Refjmf) E Fi(M'^,M) © H''{R'^,R), with the symplectic form 

a;(V', </>) = Im / 4>- 

Equation ( ll.ip is invariant under time-translation (t i— > t+to) and gauge (phase)-translation (j) 4>e''"' , 7 G 
M yielding, by Noether's theorem, the conservation laws 

Conservation of energy: <?[^/'(t)] = <S[^/^(0)]; 

Conservation of particle number (optical power): AA['(/'(i)] = M['ijj{0)], where 

Ar['0] = / IV'I^ dx. 



We make the following 

Assumptions on the potential V and nonlinearity /: 

(fA) /(r) is a smooth function satisfying /(t) = 0{t) for \x\ is small. Thus, the nonlinearity in NLS is 
cubic at small amplitudes, i.e. /(IV'P)^ ~ ^IV'PV'- 

(VA) V is smooth and decays exponentially |x| tends to oo. 

To ensure the global well-posedness of the initial value problem for (jl.ip we impose: 

(fB) Subcritical nonlinearity for large amplitudes: 

1/(01 < c(l + for some /3 G [0, f ) and 

1/(01 < c(l + |Cr~^) for some a G [0, (^^y:^)- Here, s+ = max{s,0}. 
The following well-posedness theorem can be found in [51 142] . 



Theorem 3.1. Assume that the nonlinearity f satisfies condition (fB), and the potential V satisfies (VA). 
Then equation ( \1.1]) is globally well-posed in , i.e. the Cauchy problem for Equation ( \1.1]) with a data 
Tp{0) £ has a unique solution Tp{t) in the space , which depends continuously on V'(O). Moreover, the 
solution ijj^t) satisfies conservation of energy and conservation of particle number. 
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4 Bifurcation and Lyapunov Stability of Solitons in the Weakly Non- 
linear Regime 

In this section we discuss the existence of sohtons in the weakly nonhnear regime. The following arguments 
are similar to those in \39\ [55j except that the excited states are degenerate. We assume that the linear 
operator —A + V has the following properties 



(Eigv) The linear operator —A + V has two eigenvalues eo < ei < with 2ei > cq- bq is the lowest eigenvalue 
with ground state (j)iin > 0, the eigenvalue ei is degenerate with multiplicity N and eigenfunctions 



tim clin clin 
SI ) ?2 ) ■ ■ 'i 



Remark 1. In appendix \13.1\ we construct a class of double-well examples, V, in dimension d = 3 and 
with multiplicity N = 2. 

The following result shows that nonlinear bound state solutions (i?i)^,A) of NLS/GP, (11.20 . bifurcate 
from the zero state and the linear ground state energy (0, A = — eo). 

Proposition 4.1. Suppose —A + V satisfies the conditions in {Eigv) above. Then there exists a constant 
6o > and a nonempty interval Tso C [— eo — Sq, — cq + 6o] such that for any X S Is^ has solutions 

of the form ip{x, t) = e^-^^cf)^ G with 



1 . 



= 5{\) ( + O {5{\)) ) , 6{X) = 0{\X- \eo\ 1^) (4.1) 

for I A — |eo| | small. 
Moreover, 

dx(b^{x) <ce-^l^l, (4.2) 



b^(x)\ < ce""^!^! and 



and similarly for the spatial derivatives of cj)^ and dxcfy' 



Remark 2. Suppose /(|^P)^ = —g\4>?' + o(I^P)- 

Then, for g > (repulsive case) we have forX^^^ = (— eo, — eo + So). 

For g < (attractive case), we have Isq = (— ^o — 6q, — eo). 

Finally, we conclude this section by noting that for 6' < 6q sufficiently small that soliton manifold, Ms' , 
see (|1.3p . is orbitally stable; see the discussions in the introduction and [571 [39l [26]. 

5 ^(A) = JH{X), the Linearized Operator about 0^ 

We now turn to a discussion of the operator obtained by linearization around the soliton and the existence 
of neutral modes with non-zero frequencies. Rewrite Equation ( ll.ip as 

where the nonlinear map G(^) is defined by 

G(V') = -i(-A + X + V)iP + ifM^)i^. (5.1) 
Then the linearization of Equation ( II. ip can be written as 

dG(<A^)x, (5.2) 



9x ^^,,An 



dt 
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where dG{(t)^) is the Prechet derivative of at (j)^. It is computed to be 

dG((A^)x = -i(-A + \ + V)x + ^/[(<A^)']X + i/W-fU^fix + X). (5.3) 

This is a real hnear but not complex linear operator. To convert it to a linear operator we pass from 
complex functions to real vector-functions 

V X2 

where xi = R^X and X2 = Imx- Then dG{(f>^)x < — ^ L{X)x where the operator L{X) is given by 

L(A) = JH{\) (5.4) 

where J is a skew-symmetric matrix 



J :-- 

and H{\) is a self-adjoint matrix 

H{X) := 

with 



1 
-1 



L+(A) 

L_(A) 



L_(A) :=-A + A + F-/[(0^)2] 

and 

L+(A) :=-A + X + V- f[{cf>^f] - 2f' [{<f>>^f]{ct^'f . 
We extend the operator L(A) to the complex space H'^{R'^,C) H'^{W^,C). 

5.1 The Spectrum of L(A) 

The operator L{\) has the neutral modes: 

Proposition 5.1. Let L{X), or more explicitly, L{X{d),d) denote the linearized operator about the the 
bifurcating state (/>^, A = X{S). Note that A(0) = — eo. Corresponding to the degenerate e-value, ei, of 
—A + V, the matrix operator 

L{X = -eo,5 = 0) 

has degenerate eigenvalues ±iE{—eo) = ±i(ei — cq), each of multiplicity N. For 6 > and small these 
bifurcate to (possibly degenerate) eigenvalues ±iEi{X),..., ±iEisf{X) with eigenf unctions 

±ir)i J ' \^ ±ir]2 / ' 'V ^^Vn 
with 

and 

7^ lim = lim r]j G span{^j'"', j = 1,2,- ■ -jN} in spaces for any k > 0. (5-5) 

Moreover, for S sufficiently small 2Ej{X) > X, j = 1, 2, • • •, A'', (nonlinear coupling of discrete to continuous 
spectrum at second order). 

For the case of a radial potential, V = V{\x\), the neutral modes have the following structure: 
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Proposition 5.2. // the potential is radial, V = V{\x\), then c/)^, hence d\(l)'^ , is spherically symmetric. 

clin _ £2_, 
-3 - \x\^ 



If the degenerate linear excited states are of the form = rfr^'^^dxl) for some function then 



Ej = El for any j = 1,2, ■ ■ ■, N = d and we can choose £,j,rij such that = ^Cd^^l) ^'^'^ Vj — ^^(1^1) f^''" 
some real functions ^ and rj. 

Remark 3. For d = 3, the hypothesis on the linear excited states states that these are of the proportional 
to Y^{0, (j)), m = —1,0, 1, where are the spherical harmonics of degree one. 

Proof. We sketch the proof. If V is spherically symmetric then by the uniqueness of the ground states 
and the fact — A + F is invariant under unitary transformation we have (f)^, hence d\(j)^, is spherically 
symmetric. 

We now outline a proof of the existence of and rjj with desired structure. Define a linear space 

y'= = | JGF^ J(x) = g5(|x|)}. 

By the definition, L(A) : V"^ Y^. Note that, restricted to Y^, x\) is an eigenfunction of — A + of 

multiplicity one. Application of bifurcation theory to Y"^, we prove there exists an eigenfunction {(,i,ir]i)'^ G 
Y with eigenvalue Ei. The other eigenfunctions with the same eigenvalue are obtained by noting that this 
computation can be carried out for any Xj, j = 1, . . . ,d. □ 

Based on the above discussion, we assume the following: 

(SA) Structure of the the discrete spectrum of L{X) = JH{X) 

1. ad{L{X)) consists of an eigenvalue at and complex conjugate eigenvalues at ziziE{X). 

2. The discrete subspace, corresponding to the eigenvalue is spanned by the associated eigenfunctions 

:^ ) - ( r 

3. The discrete subspace, corresponding to the eigenvalue iE{X), E{X) > 0, is N— dimensional and is 
spanned by the (complex) eigenfunctions vi, V2- ■■,vn 

4. Thus, vi, V2, ■ ■ •, are eigenfunctions which span the discrete subspace corresponding to the eigen- 
value —iE(X). 

5. Moreover we observe that Jvn are eigenfunctions of the adjoint operator L{X)* with eigenvalue 
-iE{X) 

L{X)*JVn = -JL{X)Vn = -iE{X)JVn. 

Concerning the continuous spectrum of L{X), we apply Weyl's Theorem on the stability of the essential 
spectrum for localized perturbations of J(— A) [351 EZ] to obtain 

aess{L[X)) = (— ioo,— iA] n [iX,i<x)) 

if the potential V in Equation ( II. ip decays sufficiently rapidly as tends to infinity. 

The end points of the essential spectrum are called threshold energies. 
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Definition 1. Let d > 3. A function h is called a threshold resonance function of L{X) at fi = itzA, the 
endpoint of the essential spectrum, if h ^ L^, 1^(2^)1 ^ c{x)~^'^~'^^+ and h is and solves the equation 

{L{X)-n)h = 0. 

In this paper we make the fohowing spectral assumption on the thresholds ziziX: 

Thresh^ : There exists 5', with < 5' < 5q (see Proposition I4.ip . such that for A G 2^/, 
L(A) has no threshold resonances at itzA 

In the weak amplitude limit, property (Thresh;j) can be referred to the question of whether the scalar 
operator, —A + V{x) has a threshold (zero energy) resonance. In [3^ it was shown that —A + V has a 
zero energy resonance or eigenvector if and only if the operator / + (—A + iO)~^V : (x)^L^ {x)'^L?' is 
not invertible. Moreover, this operator is generically invertible. That is, if we replace V by qV , where q is 
a real number, then we have non-invertibility for only a discrete set of q values |361 [33] 

The reduction from the properties of L{\) to those of —IS. + V is seen as follows. Let 

Then 

H := -iA*L{\)A, (5.7) 
It follows that ibiA are threshold resonances of L{\) if and only if ibA are threshold resonances of H. 
We next observe that H \s a. small perturbation of a^{—/S. + y + A), where ^3 ~ f q i ) " 1'^'^^^^, 



a computation of H yields 
where 

and for some c > 



H = Ho + Vi + VsmalU (5.8) 
//o := (-A + A)c73, Vi:=Va:, (5.9) 

\Vsmall\ < e~'l''l 0(1), 



where o(l) — > as |A — |eo|| 0. 

Therefore, the generic validity of (Thresh;j) from the generic absence of zero energy threshold reso- 
nances for —A + V by the following result, proved for d = 3 using results in [18]. The proof for general 
dimensions is similar. 

Proposition 5.3. Let d = 3. If the operator I + (—A + iO)~^V : —>■ (x)^L^ is invertible, then 

(Thresh;)^^ holds when |A — |eo|| is sufficiently small. 

Proof. We begin by proving that the operator I+{Ho±X+iO)^^[Vi+Vsmaii] ■ {x)'^L'^ (x)^L^ is invertible. 
Observe that — 2A ~ 2eo is not an eigenvalue of the operator —A + V, hence I +{—A + 2X)~^V is invertible. 
This, together with the hypothesis, implies that I+{Ho±X+iO)~^Vi is invertible with a uniformly bounded 
inverse. On the other hand the norm of the operator (Hq it A + iO)~^Vsmaii is small when |eo + A| is small. 
Hence I + {Hq ± X + iO)~^Vi + Vsmaii] = [I + {Hq ± X + iO)-^Vi][l + {1 + {Hq ± X + iO)-Vi)-Vw/] is 
invertible when |A — |eo|| is small. Moreover in [18] it is proved that the operator L{X) has no threshold 
resonance functions if the operator / + {Hq ± A + iO)~^[Vi + Vsmaii] '■ {x)'^L'^ (a;)^L^ is invertible. This 
completes the proof. □ 
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Choice of basis for degenerate subspaces 

In our analysis, it is important tliat we choose an appropriate bases of tlie degenerate eigenspaces corre- 
sponding to ziziE{X). We present this choice of basis and its construction here. 

Proposition 5.4. There exist real functions ^n, Vm n = 1,2, ■ ■ ■, N, such that 



span 



and for any m,n G {1, 2, ■ ■ •,N}, 



irin 



f U''f]{(t^^f{Urin - inr]m)dx = 0, 



span{vi,V2, ■ ■ -^vn], 



(5.10) 



(5.11) 



(0^, Cn) = {dx(t)^, r]n) = 0, {Cn, Vm) = 6m,n- 

The proof is given in the Appendix 113. 4[ 

Remark 4. is spherically symmetric, then ^„ = ^'^(l^;]) and rjn = ^r]{\x\), n E {1,2, ■ ■ ■,N = d}; 

see Lemma \5.S^) . Therefore ( \5.10\) trivially holds because (,mVn — ^n'Hm = 0. 

We conclude this section with the explicit form of the projection Pdisc, whose proof for dimension 
one can be found in [23]. The proof for general dimension is similar, and hence omitted. Recall that 

Proposition 5.5. For the non self-adjoint operator L[\) the (Riesz) projection onto the discrete spectrum 
subspace of L{\), Pmsc = Pdisc{L{\)) = P^-^^, is given by 



Pm 



disc 



5aII0^IP 



N 
n=l 







in 
iVn 





-irjr, 



+ 



dx<b^ 



in 

-illn 



(5.12) 



irjn 
in 



We define the projection onto the continuous spectral subspace of L{\) by 



P, = P,iL{X)) = P^ 



I 



Pdi 



disc- 



(5.13) 



5.2 Estimates on the Propagator 

We will need estimates of the evolution operator U{t) := e*'^'^'^^-' for Ai G X. Recall that L{\i) has two 
branches of essential spectrum [iAi,ioo) and (— zoo,— iAi]. We denote by P4. = P^^ and P- = P^^ the 
spectral projections associated with these two branches of the essential spectrum. Hence, P^^ = P+ 

Theorem 5.6. Let d > 3 and define k := [|] + 1 and u := Assume that 2E{Xi) > X, so that 

±2iE{Xi) G aess{L{Xi)). 

Then, for any time t > and Ai G T there exists a constant c such that 



(x)-'^(-A + l)l[/(t)(L(Ai) ± 2iE{Xi) - 0)-"P±/i||2 < c(l + t)-l\\{xy{-A + l)lhh 



(5.14) 
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with n = 0, 1, 2. 

For any time t £ (— oo, oo) and Xi £ I there exists a constant Cj such that 



ii(^)-(_A + i)'2U{t)p±hh < cx (1 + \t\r^\{xn-A + i)i/i||2 

\\Uit)P±h\\^<Ci \t\-l\\hh 
\\U{t)P±h\\^<Cx{l + \t\)-l{\\h\\H. + \\hh) 
\\U{t)P±hh <Cj{l + \t\)-H\\h\\H, + \\h\\i) 

\\{xy^u{t)Pihh <cj{i + \t\rimi + \\hh). 



(5.15) 
(5.16) 
(5.17) 
(5.18) 
(5.19) 



We refer the proof of the estimates to [HI |23l [52l [27|. That the constant, Cj can be taken to hold 
uniformly for Ai G T, see [T2|fTilfT3]. 

6 Matrix Fermi Golden Rule 

As highlighted in the introduction, the decay of neutral mode components, associated with the linearized 
NLS/GP equation, is necessary for asymptotic stability of the soliton manifold A4x- We shall prove that, 
after near-identity transformations, the system governing these neutral mode amplitudes is (jl.6p : 



where ziziE{X) are complex conjugate — fold degenerate neutral eigenfrequencies of L{X) = JH{X), T is 
symmetric and A is skew symmetric. It follows that 



Our strategy to show that \z{t)\ tends to zero, is based on proving that T is positive definite and that the 
corrections to (|6.2p decay sufficiently rapidly as t tends to infinity. If L{X) has a complex conjugate pair of 
simple neutral eigenvalues, then T reduces to a non-negative scalar. If L{X) has multiple, well-separated 
pairs of neutral modes, then F reduces to a diagonal matrix [371 EH [3 [Ml SSI [531 [51]. The present case of 
problem of degenerate neutral modes is more involved due to coupling among the various discrete modes 
and with the continuous spectrum. Our computation yields a non-diagonal FGR matrix, T. In this section, 
we display the expression for F, state a result on its general properties. The detailed derivation of the 
expression for F is carried out in section [lOl 

6.1 The FGR matrix, T{z,z) 

To construct F we must first introduce some notation. 
Define vector functions G^, A; = 1, 2, • • •, A^, as 




6>0 



(6.1) 



dt \z{t)\^ 



2 z* r{z,z) z + .... 



(6.2) 




(6.3) 



with the functions B{k) and D{k) defined as 



B{k) 
D{k) 



f'U^?]<P^ [ 3(z • 06 - {z ■ v)Vk ] - 2f" [(<A^)2] (0A)3 . ^)^^ 
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where 

N N 

Z-C ■= '^ZnCn, Z-rj := Znrjn- 
n=l n=l 

In terms of the column 2- vector, Gk, we define a. N x N matrix Z(z, z) as 

Z{z,z) = {Z^^'^\z,z)), l<k,l<N (6.4) 

and 

Z(''''\z,z) = -{{L{X) + 2iE{X)-0)-^P,Gi,iJGk) (6.5) 

Since Pc{L)*J = J Pc{L), a consequence of L = JH and H* = H (see Proposition 113.91 ). we have the 
more symmetric expression for Z^^'^^: 

ziKi) ^ -{{L{\) + 2iE{\)-Q)-^P,GuiJP,Gk) (6.6) 

Finally, we define T{z, z) as follows: 

V{z,-z) ■■= \[Z{z,-z) + Z-*{z,z)]. (6.7) 

Thus, 

[T{z,z)]ki = - Re {{L{X) + 2iE{X) -0)-^P,Gi,iJPcGk) (6.8) 



Concerning the properties of T, we have the following general result: 
Theorem 6.1. (Matrix nonlinear Fermi Golden Rule) 

(1) T{z,z) = r(z,z; A) is a non-negative symmetric N x N matrix, displayed in 116. 8\) . 

(2) Define 

K{X,G) := min ^i^'^} 

Then, K{X,G) depends continuously on X and G (in the space {x)'^L°°), where G = {Gi, . . . ,Gn), 
defined in i6.3\) . 



We shall require the following Fermi Golden Rule resonance condition: 

(FGR) There exists 5', with < 6' < 6 (see Proposition 14. ip and a constant C > 0, such that for any 
s = (si, • • •, snV, z = (zi,- ■■, zn)'^ G C^, 

s* r{z,z;X) s> G (6.9) 

where A G 2s'- 

Remark 5. In the weakly nonlinear regime, see section [5.11 E{X) ~ ei — cq, A ~ — cq and therefore the 
condition for resonance with the continuous spectrum at second order is: 2E{X) — A ~ 2(ei — eo) + eo = 
2ei - eo > 0. 

Our next result is a reduction of the condition (FGR), for the class of double- well potentials discussed 
in appendix 113.11 to an explicit condition on the operator V . 
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Proposition 6.2. Let V denote the double-well potential satisfying condition {Eigy), and constructed in 
Appendix \13.1[ Thus, —A + V has two negative eigenvalues, cq < ei < 0, with 2e\ — cq > 0. The excited 
state eigenvalue, 62, is degenerate of multiplicity N = 2, with spanning eigenf unctions {Ci^^i^i*"}- 
Let /(IV'P) = — Assume the matrix 

( Im({-A + V-[2e,-eo]-iO)-'Pc <t>Un e^^^Ci*", <A^in CK") (6.10) 

\ \ / / l<m,n<2 

is positive definite. Then, there exists 5' > 0, such that for the (f)^ denotes the soliton of Proposition \4-l\ 
if \ X — |eo| I < S' then K{X,G) > 0. and the Fermi Golden Rule condition (FGR), holds by taking 
C = inf^gjy X(A, G(A)) > in i6.9\) . Here, Is' denotes a sufficiently small subinterval of the range of X— 
values for which the soliton exists; see Proposition \4-l\ 

Remark 6. Positive definiteness of the matrix 16.10\) is equivalent to 

Im ((-A + V- [2ei - eo] - lOy^Pc ct>iin (^iC^ + ^2^2")', 'Pun (^iCf" + ^2?^*")') > C\z\\ 
for all zi, Z2 G C^. 

Proof of Proposition 16. 2t In what follows we sketch the proof, which is very similar to the case iV = 1, 
see imiSl]). 

Recall the transformation of L{X) in equation ()5.7p : 

(L(A) + 2iS(A) -0)-i = [iAHA* + 2iE{X) -{)]-^ 

= -iA[H + 2E{X) + i^]-^ A* 

= -iA[Ho + 2E{X) + iQ]-^A* + iA[H + 2E{X) + iQ]-^Vsmaii[Ho + 2E{X) + i^]-^A* 

(6.11) 

and 

[Ho + 2E{X) + ^0] - ^ p ^2^^^^ _ _ j . (6.12) 

On the other hand by Propositions 14.11 and 15.11 we have that in H^ space ij-^j — (j)^ — > n , — (pun and 

\\<P Wh"^ \\'Plin\\u2 

Fii"^") ^ lld'^Va ^^""'^^'"^ ^^^^ 1-^ ~ 1^0" ^ °- nonlinearity /(r) = 

r°", cr > 1 we have 

A*P,Y.^,G,=CUun\\%-' ( ;^^2.-l(^^^..^^^^.„)2 ) (1 + 0(1)) 

for some constant C G C. 

In considering KVlb . note that -A - 2E{X) ~ eo - 2(ei - eo) < and 2E{X) - A ~ 2ei - eo < 0. Thus, 
Lm{{—A + V + X + 2E{X))~^ F, F) = for any F. Furthermore, lle""^!^' ysma//||L°° is small for some c > 0, 
we have 

K{X, G) = \Cm^uJ%-^Ko{l + 0(1)) (6.13) 

with 

Ko := /m ((-A + y + eo - 2ei - iO)-' Pc{cpHnf''-\zi^f^ + Z2^'r?, (0/^n)''^-'(^id^" + ^2^2")') (l+o(l)). 

The proof is complete. In the appendix 113.31 we have a simpler formula for (FGR) when the potential V 
is spherical symmetric. 

The proof of Theorems 16.11 is deferred to Appendix 113. 2[ 
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7 Main Theorem 

In this section we state precisely the main theorem of this paper. Recall the notations ^ = (■^i, • " " t£,n) 
and r] = (r/i, . . . , t^tv) for components of the neutrally stable modes of frequencies ziziE(X) of the linearized 
operator. Recall the definition of the interval Z in (jl.4l) . 

Theorem 7.1. Assume conditions: 

(fA) (fB) on the nonlinearity /(l^/^p); section\^ 

(VA) on the potential V{x); section\^ 

(SA) on the structure of the discrete spectral subspace of the linearization about cj/^ ; section\^ 

(Thresh;^) on the absence of threshold resonances; section [57T\ 

(FGR), the nonlinear Fermi Golden Rule resonance condition, W. 9^) : section\^ 

Fix u > Q sufficiently large and let k = [|] + 1, where d > 3 denotes the spatial dimension. 
Then there exist constants c, eo > such that, if for some Aq G X 

, < c|z(°)|<eo, (7.1) 

A: ,f 

then there exist smooth functions 

X{t) : R+ ^ T, 7(t) : M+ ^ M, z{t) : M+ ^ C^, 
R{x,t) : X M+ ^ C 
such that the solution of NLS evolves in the form: 
i;{x,t) = e'^oMs)ds^n{t) 

x[ c/)^ + ai{z,z)dx^^ + ia2{z,z)cl)^ + {Re z) ■ ^ + i{Imz) ■ r] + R ], (7.2) 
where lim^^^oo A(t) = Aqo, for some Aqo G Z. 

Here, ai{z, z), a2{z, z) : x M and z — z : x C'^ — > are some polynomials of z and z, 

beginning with terms of order jzp. 

Moreover: 

(A) \z{t)\ < c(l + 1)~2 and z satisfies the initial value problem 

19 

dtz = -iE{X)z -T{z,z)z + A{z,z)z + 0{{l + t)-ir) (7.3) 

where T{z, z) is a symmetric and positive definite matrix defined in ( \6.4\ ), ^{z, z) is a skew symmetric 
matrix. 

(B) R{t) = {Re R{t) , Im R{t))^ lies in the essential spectral part of L{\{t)). Equivalently, R{-,t) 
satisfies the symplectic orthogonality conditions: 

u){R,i(j)^) = uj{R,dx(l)^) = 0, 

u{R, ir,n) = uj{R, = 0, n = 1, 2, • • •, iV, (7.4) 
where uj{X,Y) = Im J XY and satisfies the decay estimate: 

\\{l + xY''mh<c{l + t)-\ (7.5) 

Remark 7. We conclude this section by stating that all the hypothesis except (FGR) in our main result 
applies to the double-well example of appendix \13.1^ see Proposition \6.S\ for a reduction of (FGR) is to 
an explicit condition on the spectral condition on —/S. + V. We expect (FGR) to hold generically in an 
appropriate sense. 



inf 



P«7 



+ {Re z^^^)-i + i {Im 
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8 Reformulation of The Main Theorem 

In proving Theorem 17.11 we establish more detailed characterization of the perturbation about Aij. 



First, we introduce the following simplifying notation: we always use z to stand for a complex 
A^-dimensional vector z = {zi,Z2,- ■ -^zn) and an upper case letter or a Greek letter with two 
subindices, for example Qmni to stand for 



a|=m, |6|=n k=l 



N 



N 



where a,h ^ N^, \a\ := ^^a^. We refer to this kind term as {m,n) term. 



k=l 



Theorem 8.1. The following more precise decomposition of the solution in Theorem 7.1 holds. The 
perturbation R in ( \7.S^ can he decomposed as 

R= Rm,nW + R (8-1) 

m+n=2 

where Rm,n o-fe functions of the form 

Rm,n = \L{\) + iE{\){m - n) - 0]" Vm,n 

4'm,n ore polynomials of z and z with coefficients being smooth, exponentially decaying functions. The 
function R satisfies the equation 

dtR = L{X)R + M2iz,z)R + PcN2iR,z)+PcS2iz,z), (8.2) 

where 

(1) S2{z,z) = 0{\z\^) is a polynomial in z and z with X-dependent coefficients, and each coefficient can 
be written as the sum of functions of the form 

(L(A) + 2iE{X) - 0)-'=P,(/>+fc(A) or (L(A) - 2iE{X) - 0)-^'P,(/>_fc(A), (8.3) 

where /c = 0, 1,2 and the functions 0±fc(A) are smooth and decay exponentially fast at oo; 

(2) M2{z,z) is an operator defined by 

M2{z, z) := ^P,J + XP,x + X, (8.4) 
where X is a 2 x 2 matrix, satisfying the bound 

\X\ < c|z|e-^«l^l. 

(3) N2{R,z) can be separated into localized term and nonlocal term 

N2 = Lac + NonLoc (8.5) 
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where Loc consists of terms spatially localized (exponentially) function of x ^ M.'^ as a factor and 
satisfies the estimate 

||(x)"(-A + l)Loc\\2 + \\Loc\U + WLocU < c ( \z{t)f + \z{t)\ ||(x)-"(-A + l)i?||2 ) , (8.6) 



and NonLoc is given by 



NonLoc := f{Ri + R^)JR, 



(8.7) 



and consists of purely nonlinear terms in R, with no spatially localized factors. Here v is the same 
as in Theorem \7.1\ 



(4) Denote by Remainder (t) any quantity which satisfies the estimate: 

\Remainder{t)\ < |z(t)|4 + ||(^)-(_A + l)ii(t)||2 + ||i?(t)||^ + |z(t)| \\{x)-^ mh 
The functions A, 7, z have the following properties 
(A) 



(B) 



A = Remainder {t); 
7 = Ti^i + Remainder{t) 



with 



T 



1,1 



(8.9) 
(8.10) 

(8.11) 

(8.12) 

where T{z, z) is the N x N positive definite matrix defined in ( \6. 7| j, A(z, z) is skew symmetric. 

9 The Effective Equations for i, A, 7 and R 

In this section we derive equations for i, A, 7 and R. 
We decompose the solution as 



(C) the vector z satisfies the equation 

dtz = —iE{\)z — T[z, z)z + h.{z, z)z + Remainder{t) 



+ ai4>\ + ia2(t>^ + ^(a„ + Pn)in + i^i^n + qn)r]n + R 



N 



N 



J J* X(s)ds i-y{t) 



n=l 

(j)^ + ai(l)\ + ia2(p^ + {a+p)-C+{/3 + q)-V + R 



n=l 



(9.1) 



Here and going forward, we use the notations: 

a = (qi, . . . jQat)"^, P = (Pi,. . . ,(3n)'^ , 
^ = {^i,---,^nV, r/ = (??i, . . . ,7?Ar)'^. 

Introducing 

z = a + i(3, 
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we have 



1 _ 1 _ 

a = -{z + z), (5 = —{z-z), 



2' ' ' 2i 

and we seek aj = aj{z,z) = Od^l)^, j = 1,2 and p„ = Pn{z,z) = 0{\z\'^), = qn{z,z) = 0{\z\'^), 



polynomials in z and z, which are of degree higher than or equal to two. Substitution of the Ansatz ( 19.1 
into NLS, equation ( II. ip . we have the following system of equations for R with 



R 



Ri 

R2 



, Ri = ReR, R2 = ImR : 



dtR 



L{X)R + jJR-JN{R,z 

+ 
+7 



dxcP^[\ + dtai] 
(/>^[7 + (9ta2-ai] 
^.[E{X)i(3 + q)-dt{a + p)] 
-f] ■ [EiX)ia + p) + dtiP + q)] 
{f3 + q)-r] \ aidl^^ + ia+p)-dxC 

-ia+p)-^ J [ a2dx(l)^ + if3 + q) ■ dxv 



(9.2) 



Here, 



N 



Re N{R,z) 
Im N{R,z) 



, JN{R,z) 



Im N{R,z) 
-Re N{R,z) 



(9.3) 



ImN{R, z) := /(|0^ + h + ih\^)l2 - m^f]h 

ReN{R,z) := [/(|0V /i + ^/sl') - f[{<P^fW + h) - 2f' [{cl,'f]{ct>'fh 
h = A1+A2 + R1, h = Bi + B2 + R2 
Ai = a-i, A2 = aidx<p^+p-i, 
Bi = P ■ r], B2 = a2(t)^ + q-T). 

Prom the system of equations ( I9.2p and the orthogonality conditions ( 15. lip we obtain equations for 
A, 7 and = On + ifin, n = 1, . . . ,N : 



dt{an+Pn) - E{X){f3n + qn) + {ImN {R, z) , r]n) = Fin, 

dt{Pn + qn) + E{X){an+Pn)-{ReN{R,z),^n)=F2n; 



7 + dta2 - ai 
X + dtai + 



1 



{^\dx^^) 
1 



{ReN{R,z),dx 



{ImN{R,z),(t>' 



Fa. 



{(p\dx<P^) 

Finally, the scalar functions Fj n, j = 1,2, F3, F4, are defined as 

Fin = 7((/5 + q) ■ V,Vn) - Xai{dl(f)^,r]n) - X{{a + p) ■ dxC,r]„) -^{R2,r]n) + X{Ri,dxr]n), 

F2n = -7((a +P) • (,^n) - Xa2{(t>x,U) - HiP + q) ■ dxri,in)+i{Rl,in) + KR2,dxin), 

F3 = 



(9.4) 
(9.5) 
(9.6) 

(9.7) 



1 



1 



X{R2, (AL) - 7(^1, 0a) - (7(a + P) • e + Aa2(AA + A(/3 + q) ■ dxv, 4>\) 
X{Ri,<pl) + 7(i?2, (A^) + (7(/5 + q)-r]- Xaidl^^ - X{a + p) ■ dx^, </)^ ) 



(9.8) 
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Remark 8. (a) Recall the estimate of Remainder in ( \8. <gj) . By Equations ( \9.4^ -( \9. 7| j we have 

A, 7, dtZn + iE{\)zn = 0{\z\'^) + Remainder. (9-9) 

(h) The functions aj{z,^),j = 1,2, qn{z,z),n = 1,...,N will be chosen to eliminate "non- 

resonant" terms: 2 < \m\ + |n| < 3. 

Finally, we derive an equation for 

R = P^(t)R = P^R, 

the continuous spectral part of the solution, relative to the operator, L{X{t)). Applying Pc = Pc^^^ to ( I9.2p 
and using the commutator identity: 

PcdtR = dtR - XdxPcR 

we obtain 

dtR = L{X{t))R - P^^'hN{R,z) + L^^^^^R + Q. 



The operator L^j^ and the vector function Q are defined as 



L 



= xidxP,^^'^)+jP,^^'h, 

p\(t) ( i{(3 + q)-r]- Xaidl 



A(a + p) ■ 5aC 



(9.10) 
(9.11) 



- X{P + q) ■ dxv 



(9.12) 
(9.13) 



We now summarize the preceding calculation in the following: 



Proposition 9.1. (Reformulation of NLS) 

Using the Ansatz ^9. 1]) 



J Jo H^)ds„i-yit) 



+ ai0^ + M2</)-^ + {a+p)-£, + {P + q)-r] + R 



(9.14) 



NLS can be equivalently expressed as a coupled system of equations ( |9.^P^?.5P9.6]P. 7\ ) for modulating solitary 
wave parameters: X{t), ^{t), neutral mode amplitudes Zn{t) = a„(t) + i/3„(t), n = 1, . . . , N and an equation 
\9.11\) , governing "dispersive part", R{t), which evolves in the continuous spectral subspace of L{X{t)) , i.e. 
P^^^^R{t) = R{t); see [SJE) . 

Finally, the functions aj = aj{z,z),j = 1,2, p{z,z),q{z,z) = {pn,qn)n=i,...,N are 0{\z\'^) and are polyno- 
mials chosen (in what follows) to eliminate "non-resonant" terms of the form z°'z'', 2 < |a| + |6| < 3. 

9.1 Extracting the 0{\z\'^) part of R{t); Proof of Equation ( lOl) 

For fixed z{t) G C^, the equation for R{t) is forced by terms of order ©([^(t)^); linear terms are removed 
due to the equations satisfied by z{t) = a{t) + if3{t). In our analysis, we need to explicitly extract the 
quadratic in z,z part of R{t). 

Thus, we consider the quadratic terms generated by the nonlinearity: 

^ ^ m,n 
ATRe 

JNm,n = JiV2,o + JiVi,i + JiVo,2 

fYi -[■- Yl — 2 

/ 2f'[{4>^f]4>^A,Bi 



m,n 



3/'[(0")']0" + 2/"[(c/>^)2](0^^3 



f[{cf>^f]cf>^Bf 



where Ai = a ■ Bi = (3 ■ rj. 
Now we have the main theorem 
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Theorem 9.2. Define 

Rm,n := [L{X) + iE{X){m - n) - 0]-^P,JNm,n, (9.15) 

and decompose R{t) as 

R = + ^ (9.16) 

m+n=2 

The function R{x,t) satisfies ( \8. S^) . 
Proof. We observe that to get 

dtR = L{\)R + L^^^.^R+ ^(A,7)^-,„ + a 



m+n=2 

- Y [dtRm,n + iE{X){m-n)Rm,n]- PcJNy2 
m+n=2 



(9.17) 



where we used the definitions of Rm,n in ( I9.15|) . the operator L^j^ and the term Q are defined in ( 19. lip . 

JN>2:=JN{R,Z)- Y J^rn,n. 

m+n=2 

Next we further decompose JA^>2 and find M2, S2 and in ( 18. 2p . We consider the functions 
JNm,n, m + n = 3, the third order terms of JA^>2 : 

E M.. ^ ^ ( E ^ H ^: ) ) + ( "^^^ ) (^'i^) 

m+n=3 \ m+n=2 ^ ' / ^ x ^ ^ / 



where, recall the definitions of Bi, A2, B2 from (j9.3p . 

GMlBl) := f'U^)\Al+Bl) + 2f"[{4>^)\4>^fAl 

and 

G2{AlBl) := [f'U>^f] + 2f"[{<t>^f][ct>^f]{Al + Bl) 

4 
3 

and X is a 2 X 2 matrix of order Izl, defined as 



+ [2/"[(,^^)2](</>A)2 + 4/"[(^A)2](^A)4]^2 



X - r ^ X - I 2/'[(0')']0'i?i 2/ [(<^^)2]<^AAi \ 

We define the linear operator M2{z, z) as M2(z, z) := X + L^\^ which satisfies the equation ( [8 
The function ^2 , in the statement of Theorem 17. H is defined as 

S2{z,z) := Y ^(A,7)^rn,n + ^ - ^ [dtRm,n + iE{X){m - n)Rm,r. 
m+n=2 m+n=2 

+ Y J^rn,n- 
m+n=3 

By Equation ( 19. 9p and 

[dt,{L{X)±iE{X)-d)-^] 
= dt{L{\) ± iE{\) - 0)-i - (L(A) ± iE{\) - Q)-^dt 
= A(L(A) ± iE{\) - {))-^dx[L{X) ± i^;(A)](L(A) ± iS(A) - 0)-^ 



(9.20) 
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we have that S2{z, z) satisfies the estimate in part 1 of Theorem 18.11 For the details we refer to [23]. 
Lastly, we define the nonlinear term 



N2{R,z) :-- 



^ JN{R,z) 



m+n=2,3 



(9.21) 



Using the smoothness of the nonlinearity /[•] and the removal of Od^p) and 0(|2;p) terms, we have that 
N2{R,z) = Loc + NonLoc, see ()8.5p . satisfying (j8.6p and (|8.7p . The computation is straightforward, but 
tedious, and is therefore omitted. 

Collecting the various definitions and estimates above we have ( 18. 2p . □ 



9.2 z(t) dependence of equations for \{t) and ^{t) 



In this section we present the proofs of equations ( 18. 9p and ( I8.10p , crucial to controlling the large time 
behavior. 

Here's the idea. Central to our claim about the large time dynamics of NLS, is that the solution settles 
into an asymptotic solitary wave, 0^°°, where \{t) — > Aqo- We show this by establishing the integrability 
and uniform smallness of A. Since we expect the neutral mode amplitudes, z{t), to decay with a rate t~2 ^ we 
require that there be no 0{\z{t)\'^) in the equation (|9.7p : A(t) + dtai{z,z) = .... The strategy is to choose 
the quadratic part of the polynomial ai{z,'z) so as to eliminate all quadratic terms non-resonant terms. 
The latter are terms whose z-behavior is: (-Zfe)^ or (^)^ and are oscillatory with frequencies ~ ±2iE{X). 
But what of the terms of the form: ZkZ^, which are resonant (non-oscillatory)? This is where we use the 
choice of basis for the degenerate subspace; section 113.41 A consequence of this choice is that there are 
no resonant quadratic terms appearing in the equation for A! The calculation is carried out below; see 
Lemma 19.41 

In what follows we use the notations N^^, ^mn stand for functions satisfying 



JNr, 



We define the polynomials ai, a2, Pk and qu, k 

ak{z,z) : 

Pk{z,z) ■■ 



1,2,- ■■,N in dSl]), see also (EZ]), as 

^m,n(-^)) k = 1,2, 



m+n=2,3,m^n 

^mliW^ k = 1,2 
m+n=2,3 



with the explicit forms 



qk[z,z) 



2zE(A)45 



1,2, 



(9.22) 



m+n=2,3 



where, recall Ti i from (j8.1ip . 

-2iii;(A)4')+4;) 

-i^(A)4'J+ylg 



1 



(9.23) 



1 



KiV27,<A^)-§Ti,i(z.7?,c/>^)] 



1 



(9.24) 
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and 



^t} '■= A^fc forn = 1, 2, k + l = 2,3, k ^ I; 



-2i£;(A)g(^)+i^(A)pg 



-3iii;(A)Q(;;)+P(A)Pi"o^ 



(9.25) 



AT 



fc=l 



k=l 



(9.26) 



p(") p(") /-)(") /-)(") 

-^fc,/ ' '^Lk ■~ '^k,l ■ 



The following is the main result. 



Proposition 9.3. Define the polynomials ai{z,z), a2{z,z), pn{z,z), qn{z,z) as above. Then, ( W. 
( \8.10\) hold and moreover 

dtX = Remainder {t), 

dt-f = Ti,i + Remainder{t) (9.27) 



dtZn + iE{X)Zn = - ( JA^2,1 

N 



+ Iti,! E / f ^ J" ) , f J ) ) + Remainder{t), (9.28) 

m=l 

where Ti^i is defined in i8.11\) and there are no zz terms in Remainder (t). Moreover, 

\Remainder{t)\ < \z{tr + \\{x)-'^{-A + l)R{t)\\l + \\R{t)\\l + \z{t)\ \\{x)-'^R{t)h (9-29) 
Before proving the proposition we state the following key observation. 
Lemma 9.4. 

(iV/-,</>^)=0. (9.30) 

Proof. Recall that 

Ai = a- ^ = -{z ■ ^ + z ■ Bi = l3-r] = —{z-r]-z-ri) 
The explicit form of JN2,o + -''^^1,1 + JNo,2 in ( 19.15^ imphes that 
Ni;^ + + Ni:^ = 2f'[{cP^f]cP^AiBi 

N N N N 

n=l n=l m=l m=l 
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By taking the relevant terms we have 

Af TV N N 

n=l m=l n=l m=l 

N N 
n=l m=l 

which together with ( IS.lOp yields 

N N 
n=l m=l 

This completes the proof of Lemma 19.41 □ 
Proof of Proposition 19.31 

Recall the estimate of Remainder in ( 18. Sp . We put ( 19. 6p ( 19. 7p in the matrix form 

\Id + M{z,R,p,q)\{ ^ ^ ]= VI + Remainder (9.31) 
V 7 - / 



with the matrix Q. defined as 



^^x;^[(i?eA^, 5a 0^) - ^Ti,! ( (z + z) • dx(t^^ )] - Ti,i - dta2 + ai 



the term Remainder is produced by ^^a^q^^^^a^ ^ ^^^^A^^^ ^{'^c^^)^'^^ ^ ) , Id is the 2 x 2 identity matrix, 
M(2;, g) is a vector depending on z, R,p and g and satisfying the estimate 

\\M{z,R,p,q)\\ = 0{\z\) + Remainder. (9.33) 

Now by the definitions of ai and 02 in ( I9.22p . we remove the lower order terms in z, z from {ImN, 
|Ti,i {{z-z)-r,,^^) and {ReN, dx(t>^) + ^Ti,! ( (z + z) • dx<t>^) to get 

17 = Z)i + 

with 



where we used ( I9.30p ; 



Now we claim that 



m+n=2,3 

{ReN- Y N^%^dxcP^) 

\ m+n=2,3 j 



dtJ^^^n^iE{\){m-n)A^^^n 
1(2) , „^ .(2) 



m+n=2,3 V dtJ^rn,n ^ iE{X){m - n)J^m^n J 

D\^D2 = Remainder. (9.34) 
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If the claim holds then estimates ( 18. 9p and ( IS.lOp follow from the estimates ( I9.3ip , ( I9.33P and the facts 
n = Di + D2. 

In the next we prove claim (j9.34p together with (|9.28p . 

By the fact we removed all the second and third order terms of JN we obtain Di = Remainder. Recall 
the estimate of Remainder in (jS.Sp . 

To estimate D2 we have to start with studying the equation for z. By the fact that 

dtZn + iE{\)zn = Od^l^) + Remainder 

in ( 19. 9p we obtain D2 = Od^p) + Remainder, hence 



A = 0{\z\^) + Remainder 

7 — Ti_i = 0{\z\^'') + Remainder 



(9.35) 



which together with the expansion of JN in ( I9.2ip yields 



dt{an+Pn)- E{X){l3n + qn)+ {N^'J' , i]n) = - -T {z - z) ■ T], r]r, )+ Remainder ; (9.36) 

k+l=2,3 

dtiPn + qn) + E{X){an+Pn)- Yl i^^t , ^n) = -^T {z + z) ■ ) + Remainder (9.37) 

k+l=2,3 

where, recall the real function Ti 1 from ( l8.1Up . Choose p„ and g„ as in ( I9.25P to remove the lower order 
terms as in the equations of A and 7, which together with the definition Zn = ctn + if3n enables us to obtain 



iz ■ T] \ 1 r]. 



^tZn + ^E{X)z^ = -(^JN2, + \T,,y^_^ ;'V-.e. 

+ Dj,{n) + Remainder (9.38) 

with D'i{n) defined as 

D3{n):=- idtPS+iik-l)E{X)Pj^^]-i ^ ^^[1 + ^ik - mX)Qt']]. 

fc+Z=2,3 k+l=2,3 

We claim that this together with the equations for A in ( I9.3ip implies that 

\D2\, \D3{n)\ = Remainder. (9.39) 

Indeed, by (j9.9p we have dtZn + iE{X)zn = 0{\z\'^) + Remainder which together with the equation for A 
in (|9.35p implies = 0(|zp) + Remainder, in turn we have an improved equation for Zn as 

dtZn = —iE{X)zn + 0(1^1^) + Remainder. 

Using this and repeating the analysis we find there is no Od^l^) term in D2 and D3, hence ( I9.39P holds 
which leads to ( I9.28P and (|9.34p . The proof is complete. 

□ 
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10 Proof of the Normal Form Equation ( 18.121) 



Recall the definitions of the functions B{n) and D{n) after ( I6.3p . Then the function JN2fl in ( I9.15P 
admits the form 

f Bin) 

n=l 



D{n) 



(10.1) 



The following is the result establishing the desired normal form of the differential equation for the 
neutral mode amplitudes, z{t). 

Theorem 10.1. By defining the polynomials ai,a2, Pn, Qn, n= 1,2, ■ ■ ■, N, as in ( \9.22\) -f \9. 25\) . Equation 
( [BAWi holds. 

Proof. Recall the definitions of JN^^m m + n = 3, from ( I9.18P and the equations for Zn in ( I9.28P whose 
first two terms on the right hand side admit the expansion 



with 



Ki{n) := - ( Xo,ii?2,o, 



/ N \ 



K2{n) 



X 



1,0 



k=l 

N 



V fc=i 



k=l 



-^2,0,-^^0,1 



Vn 



N 



\ 



k=l 

N 



\ k=l 



J 



where, recall the definition of X in ( I9.19P and we divide it into two terms X = Xi^ + Xq^i 



and Xo,i := Xi^; 



z-i r\ -ie 



+( wm']+^f"m'm']{\z-e + \z-vn '"'l' , 



2 I z-r] 




-z< 



Tin 



i^„ 



K,{n) := - Ti,i 



K^{n) := - 

where, recall the real function Ti^i from ( I8.10p 



— iz ■ r] 

z< 



^1,1; -^^1,0 



rin 
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Next we study Kj{n), j = 1, 2, 3, 4, 5. 

We start with the most important term Ki{n). Recall the definition of G„ in ( 16. 3|) . By direct 
computation we obtain 

^5,.(\)=-'^(;^{:j) (10^2) 

which together with ( I9.15P and ( 110. ip implies that 

N 

Ki{n) = Zk{{L{X) + 2iE{\) - Oy'PcGk, iJG„). 

k=l 

Define Z{k,n) := -{{L{X) + 2iE{X) - Q)-^PcGn,iJGk) and a iV x TV matrix 

T{z,z) ■.= [A{k,l)]w-\th A{k,l) ■.= ]^[Z{k,l) + Z{l,k)], l<k, I < d. (10.3) 

5 

For Kj{n) we claim that it can decomposed into the matrix form 



i=2 



Kj{n) = {Sin, 1), S{n, 2), • • •, S{n, d))z with S{k, I) + S{1, k) = 0. (10.4) 

i=2 

Define a, N x N skew symmetric matrix 

Aiz,z):= [A(i, k)] with A(j, k) := Sk,i + \[Z{k,l) - Z{l,k)]- 

This together with ( I9.28P and ( 110. 3p yields the equation for z in ( I8.12p 

What is left is to prove ( 110. 4p . To avoid tedious, but simple, computations, we only analyze part of 
K2{n) and i^3(n). 

(A) The part of K2{n) we study is 



the analysis of the other terms are similar. By ( 110. 2p we rewrite ^2,i{n) as 

4>A.^" \ / D{n) 



= -4iA^^l{dx^\ D{n)) + 4i4'J(</.^ B{n)). 



Equation i W?2^ relates {N^fl,(t)\) and {N^'^^cf)^) to ^2^0 and which together with the expression 
of JN2fl in Equation ( 110. ip yields 



N 



k=l 
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with 

{B{k),cl>^){<p\B{n)). 

By straightforward computation we have 



^{n,k) + ^{k,n) = 0. (10.6) 

By ( 110. Sp and ( 110. 6p we complete the proof for "i>2,i{'n). 
(B) To simphfy the notation we introduce p and lo 

N 



n=l 
N 

U := = ^'^ZnT]n. 

n=l 

This imphes that 

TV , TV 



n=l n=l 



TV , N 



N 



k=l 

with 



Hn,k) := i([/'[(0A)2] + 2/"[(0Y](0')'](p^fc-^%),(/^en-u;^n)) 
+^([3/"[(</'")'](0')' + 2/"'[(0^)2](</.^)4]pefe,pU 



Immediately we have 



(10.7) 



P^ = \Y^ ^ninP, ^'^ = 
n=l 

By the definition of K^{n) it is not hard to get 

^3(n) = ^^fc$(n,A;) (10.^ 



^>(n,/c) + $(A;,n) = 0. (10.9) 
This together with ( llO.Sp completes the proof for Kz (n) . 

□ 



11 Proof of the Main Theorem 17.11 

For simplicity, we present the proof of Theorem 17. II for the case d = 3; the proof can be easily modified to 
cover d > 3. The main difference is that, in controlling ||-R(t)||Loo(Kd) by ||-R(i)llH'={M'^)> for (i = 3 we take 
k = 2, while in general we need A; = [^] + 1; see section 15.21 
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11.1 Estimation strategy 

In this subsection, discuss our strategy for studying the large time behavior of solutions. 

We begin by introducing a family of space-time norms, Z(T), TZj{T), for measuring the decay of the 
z{t) and R{t) for < t < T, with T arbitrary. We then prove that this family of norms satisfy a set of 
coupled inequalities, from which we can infer the desired large time asymptotic behavior. Define 

To:=|zW|-\ (11.1) 

where z^^^ and the constant u are defined in Theorem 17. li 



Family of Norms: 



Z{T) ■, 


= max(To + t)2|z(t)|, 


7^l(T) 


7^2(^) 


■.= m^{To + t)\\Rit)\\^, 




7^4 (T) 




7^5(T) 



■■max{To + t)\\{x)-''R{t)\\H2, 
■ma^{To + t)^\{x)-'R{t)\\2 



(11.2) 



Remark on choice of norms: It is clear that a combination of H^, spatially weighted and L°° 
norms of R{t), as well as a bound on |2(t)|, are plausible choices of norms to control the large time be- 
havior. This accounts for the definitions: Z{T),TZi,TZ2(T) and 7^4(T). Our list of norms also includes 
estimation of the time decay of ||-R(f)||3 (T^s) and the local norm of an auxiliary function, R{t) (T^s). 
Why these two additional norms? As will be seen, the C{t) = \z(t)\ satisfies an equation of the form 
^ ~ —K^^"^ + c(t), where c{t) consists of various coupling terms (products) involving neutral mode ampli- 
tudes {z{t)), the ground state (cj)"^^^^) and dispersive terms {R{t)). First, neglecting c(i), we observe that 
C(t) ~ t~2. To treat c(t) as a small perturbation for large t, it is necesssary that it decay more rapidly 
than the term C^{t) ~ t~2. Without any further decomposition of R{t), we find among the coupling 
terms one of order \z{t)\ R{t)\\2. The expected decay rates of each factor, imply this term is of 

order t~2 for large t, which is of the same order as C^(^)- The resolution is to expand R{t) as a leading 
order part consisting of terms Rm.n = z'^z'^, m + n = 2 plus a more rapidly decaying correction R{t), 
with II (x)"'' -^(*)l|2 = 0{t~^~^), 6 > 0; see equation (|8.1[) . This modification yields an equation with an 
improved correction term of order |2:(t)| ||(x)~'^ -^(*)l|2 = t~2~^ ^ 6 > 0, which can be treated as a small 
peturbation in the large time dynamics. 



Remark on the estimation strategy; see also, [71 I48| . Estimation of lZj{T) proceeds by using 
the DuHamel, equivalent integral equation, formulation relative to a fixed operator linearized operator 
-L(Ai), where Ai = A(T) and T > is fixed and arbitrary. Namely, 

L{Xi)P^'R+ P^' {L{X{t))-L{Xi))R + ... 

= e^(^^)*i?(0) + r e^(^^)(*-^) (• • •) ds (11.3) 
Jo 

We can therefore apply the time-decay estimates of Proposition 15.61 to obtain bounds on local decay and 
L~ norms of Pc'Rit). However, we need bounds on R{t) = Pc''^^R{t). Since R{t) = P^'Rit) + PdL^{t), 



dtR = L{Xi)R + .. 

=^ dtPtR = 
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it suffices to bound P^lg^Rit)- Tliis is done as follows. 

3A1 _ lDHi)\r>U\ fU^r.„.,c^ p^W 



(P^sc - PZDPit)^ (because PZ^Rit) = 0) 

c-p^bpkcPit) + (pkc-p^SbptPit) 

I - (P^c -Pi!sl)) P^cR = (P^c - P^l) Pc'Rit) (11-4) 



Therefore, 

PdLPit) = (I- Ai))-i 6{\, Ai) P,^-R{t) (11.5) 
and we estimate R{t) in either a local energy, ff^(M'^; {x)^^dx) or L°^(M'^) via 

< \\PtR{t)\\x + \\PtR{t)\\x. (11.6) 

Here, 5(A, Ai) = -Pj^^^ — P^jfl is finite rank and of small norm proportional to jj^ |A(,s)| ds. 

We now derive the integral equation for P^^ , which is the basis for our time-decay estimates. If we 
write L{\{t)) = L(Ai) + L{\{t)) - -L(Ai), then Equation i WU\ for R takes the form 

dtP^R = L{\i)P^^R+ (A - Ai + j)Pc'JR + ■■■ 

Recall that L{X) has two branches of essential spectrum [iX,ioo) and (— ioo,— iA], we use P+ and P_ to 
denote the projection operators onto these two branches of the essential spectrum of L{Xi). Then we have 

Lemma 11.1. For any function h and any large constant > we have 

{x^i-A + 1) ( P^^J- i{P+ -P-))h < c II {x^i-A + l)h |L . (11.7) 



For d = 1 the proof of this lemma can be found in [7J, the proof of d > 3 is similar, hence omitted here. 



Equation ( 19. lip can be rewritten as 

dtP.^^R = L{X,)P,^-R+[^ + X-X,]iiP+-P-)R .-^^g^ 

+P^^OiR + P^^Pc^^G -P^^Pc^^hN{R,z), ^ ■ ^ 

where Oi is the operator defined by 

Oi := AP,A + ^(A) - L(Ai) + ^P^J - [7 + A - Ai]i(P+ - P_). (11.9) 

By Equation ( Ill.Sp and the observation that the operators P+, P- and L{Xi) commute with each other, 
we have 

+ ^{t~s)L{X,)+a{t,s){P+-P.)pX, p^jl ^ pXg _ pX J ^ [R^ z)]ds, ^ ^ ^ 

with a(t, s) = i [7(1") + A(r) — Ai] dr. We observe that -P+P- = P-P+ = and for any times ti < t2 
the operator 

ga(t2,ii)(P+-P-) = ea{t2,tl) p^ + g-a(t2,ii) p_ ■ R'^ ^ 

is uniformly bounded. 



We conclude this subsection by recording the following result which is used repeatedly in our estimates: 
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Proposition 11.2. Let Tq > 2. There exists a constant c > 0, such that 

* 1 1 , c ,3, , , 

: ;t rT—-^ ds < a e [0, -], (11.11) 

{l + t-s)2 [To + sf [To + tf 2 

/ {t- s)-^To + s)-^ds < c (ro + t)-hog(ro + t). (11.12) 
Jo 

Similar versions can be found in many literature, for example [47^ [7]. The proof is given in the 
Appendix 113.61 

11.2 Estimate for 7^l(^) := max{To + t)\\{xy'' R{t)\\Hk 
Proposition 11.3. 

ni<c(^To\\{xyR{o)\\H2 +nln2 + + Tp[z'^ + zni + UiUl] ^ . (11.13) 

With a view toward proving the time decay estimate of Proposition 111.3^ we now first give appropriate 
norm-estimates of the latter terms in equation (jll.Sp . 

First, the norm definitions ( lll.2p and Lemma lll.ll we estimate the OiR and Q terms 

\\{xY{-^ + i)0iR\\2 < c(^o + ^)-iz7^l, (ii.u) 

\\{xr{-A + l)g\\2 < c(To + t)-iz3. (11.15) 
Next, we estimate the nonlinear term, JN, via 
Lemma 11.4. 

||(-A + l)JNiR, z)\\, + ||(-A + l)JNiR,z)\\2 

< c(To + t)-^[nln2 + z^] + c{To + t)-l[zni +71^711]. 

Proof. Recall the definition N2{R, z) := —JN{R, z) + X]m+n=2 3 J^m,ni in ( I9.2ip and the decomposition 
N2 as the sum of Lac and NonLoc in ( I8.5p . By the fact JNm,n-, m + n = 2, 3, are localized functions we 
have the estimate 

IK- A + l)[JN{R,z) - NonLoc]\\i + ||(-A + l)[JN{R,z) - NonLoc]\\2 



< c\z\{ \z\ + \\(x)-''R\ 



< c 



(^o + ^)-lz2 + (^o + ^)-2Z7^l 



More challenging is the term NonLoc, defined in ( 18. 7p . which is purely nonlinear, having no spatially 
localized factors. We use the estimate 

||(-A + l)NonLoc\\i + ||(-A + l)NonLoc\\2 

^ C Q|fi||^2 ||-R||oo + ll-Rlli^z 11-^11^ ) 

< c{To + ^)-l7^|7^2 + c{To + ty^TZiTZl 
by the fact /(x^)x is of the order around x = for d = 3. 

Collecting the estimates above we have ( I11.16P . □ 
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Proof of Proposition 111.31 By Equation ( Ill.lOp . Estimates ( I5.15P and ( I5.19P for d = 3 we have 



||(x)-(-A + l)P,^^i?(t)||2 

< ||(x)-^(-A + l)e*^(^i)P,^ii?(0)||2 

+ 11 lii^yi-^ + l)e(*-^)^(^i)P,^i [Oi{s)R + P^g - P^JN{R,z)]ds\\2 

< c(l + t)-i||(x)-(-A + l)i?(0)||2 

+/o*(i + 1 - s)-l \\{xy{-A + i)[OiR + p^g]ds\\2 

+ /o*(l + t - s)-i(||(-A + l)P,^JNiR{s), z)h + ||(-A + l)P,^JNiR{s), z)\\2)ds. 
Therefore by estimates ( lll.l4p -( I11.15P and ( I11.16P we have 
\\{x)-^{-A + l)P,'-R\\2 

< c[ii+tri\\{xn-A+i)Rm\2 

+ /o (1 + t - s)-f (To + s)-^ds 7^l7^2 + Z^ + TP [Z^ + ZTli + 7^47^2] 

Using the time convolution estimate (|ll.lip we obtain 

\\{x)-^{-A + l)PtR\\2 
< c{To+t)-^ 

This implies Proposition 111.31 



(11.17) 



ro||(x)^i?(o)||2 + 7^|7^2 + z^ + t^^ ^[z^ + ztzi + n^nl] 



□ 



11.3 Estimate for 7^2(^) := max{To + t)\\R{t)\\, 



Proposition 11.5. 



7^2 < c 



To ||i?(o)||i + ||i?(o)||H2 + z2 + 7^27^2 + To ^[z3 + z7^l+7^?] 



(11.18) 



To prove this, we use the following result, whose proof very similar to that of Lemma 111.41 ^ind hence 
omitted. 

Lemma 11.6. 

\\P,^JN{R,z)\\^ + \\P,^JN{R,z)\\j,2 

< c(To + ^)-l[z2 + 7^27^2] + c(^o + ^)-i[z3 + z7^l + 7^2]. 

Proof of Proposition 111.51 By Estimate ( I5.17P for (i = 3 and Equation ( Ill.Sp we have that 

< ||e*^(^i)Pc^ii?(0)||oo + /o \\e^^~'^^^^'^ P^'[Oi{s)R + P^g - P^JN{R,z)]\\oods 

< c(l + t)-t(||i?(0)||i + ||i?(0)|U.) (11.20) 
+c/q*(1 + t - s)-i [\\0,{s)R + + \\0,{s)R + P.^gMds 

+cJ^il + t-s)--2{\\P,^JN{R,z)\\i + \\P,^JN{R,z)\\H2)ds. 

By the properties of Oi (Equation ( I11.9P ) and g (Equation ( 19. lip ) we have 

\\Oi{s)R + P,^g\\i + \\Oi{s)R + P,^g\\H2 < c{To + t)-t + z% 
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This together with ( Ill.lOp yields 

WP^^mWoo < c{To + t)-'[To\\R{0)h+To\\R{0)\\H2+Z^ + nln2 

+tP[z^ + zui+nl]]. 

This imphes the proposition. 



11.4 Estimate for U^iT) := max + \\R{t)h 



□ 



Proposition 11.7. 



7^5 < c 



To ||i2(o) II i + ||i?(o) 11^,2 )+z' + To~'[nin2 + z' + zni + ni + ni 



(11.21) 



We use the following lemma to prepare for the proof 
Lemma 11.8. 

3 



||jA^(i?,z)||3 <c{To + t)-^z^ + c{To + t)-2 [nln2 + z'^ + zni + nl + nl ] . (11.22) 

Proof. As in the proof of Lemma 111.41 we decompose JN into the localized term JN{R, z) — NonLoc 
and non-localized term NonLoc. The estimate of the first term is similar to that of Lemma 111.41 hence 
omitted. The nonlocal term NonLoc defined in ( 18. 7p admits the estimate 



NonLoc\\a < c I / \R\i j < c ||i?||i||i? 



I oo ■ 



By using the definitions of estimating functions on all the terms above we have ( lll.22p . □ 
Proof of Proposition 111.71 By Estimate ( I5.18P for d = 3 and Equation ( I11.16P we have that 

WPc^^mh 



< ||e*^(^i)P,^ii?(0)||3 + /o ||e(*-^)-^(^i)P<,^i[Oi(s)i? + P,^g - P^JN{R,z)]\\3ds 

< c(i+r^(iii?(o)iii + iii?(o)iiH2) 

+c /o - s)-^ ||Oi(s)^ + P^QUds + /g (t - s)-3 \\P^JN{R, z)\\-^ds. 



(11.23) 



By the properties of Oi (Equation ( I11.9P ) and Q (Equation ( I9.1ip ) we have 

\\Oiis)R + P^g\\ 3 < c(To + t)-l [ZUi + Z\ 

2 

This together with ( I11.19|) and ( I11.12P implies 

||^'^ii(t)||3 <c(ro + t)-hog(To + t)[To^||i?(0)||i + ro||i?(0)||H2 

+ z^ + Tq"^ [7^^7^2 + + zTZi + nl + nl]]. (11.24) 

This estimate and the definition of TZ^ yield the proposition. 

□ 
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11.5 Estimate for 7^3(^) := max(To + t)5 ||(x)-^^(t)||2 



t<T 



Proposition 11.9. Let the constant v the same as in ( \5.14\ )-( \5.15\) with d = 3. Then 



TZs < c 



3 



{xrRm2+ 



(11.25) 



As usual we estimate the nonlinear term N2{R, z). 
Lemma 11.10. 



/o 1 1 (x) -'^e(*-^)^ W P,^iV2 (i? , ^) 1 1 2 



< c (To + t)-l Tq 2° [ z3 + zui + 7^3 + 7^27^4 ] . 



(11.26) 



Proof. We start with the function Recall that = Lac + NonLoc in ( 18. 5p and the estimate of Loc 
after that. 

The nonlocal term NonLoc defined in ( 18. 7p admits the estimate 



\\NonLoc\\i + \\NonLoc\\2 < c 



+ \\R\\r 



< c 



l-^lli + ll-^llLll-RIb 



By the definition of estimating function we have 

\\NonLoc\\i + \\NonLoc\\2 < c{To + t)-l [login + t) TZl + {To + ty^TZlTZ^ 

< ciTo + t)-lT^^' [7^i + 7^i7^4 ] . 

Finally we prove Equation ( I11.26P . By the propagator estimates ( I5.17P and ( I5.19P we have 

/o \\{x)-'e('-'^^Wp,^^P,^N2iR, z)\\2ds 
< c/(J(l + t-s)-i \\\NonLoc{R,z)\\i + \\NonLoc{R,z)\\2 + \\{xyLoc\\2 

This together with the estimates of Loc and NonLoc above yields ( lll.26p . 



ds. 



□ 



(11.27) 



Proof of Proposition 111.91 By the same techniques as in deriving Equation ( I11.8P we have the 
following equation 

dtP^^R = L(Ai)P,^ii?+(7 + A-Ai)i(P+-P_)fl 

+P{z, z)R + PtS2{z, z) + P^^P^N2{R, z), 

where the operator P{z, z) is defined as 

P{z, z) := PtM2{,z, z) - (7 + A - Ai)i(P+ - P_) + P^^ (L(A) - L(Ai)), 

and the terms £'2(2,2), M2(z,z) and P^N2{R, z) are defined in Theorem 18.11 

Rewrite Equation ( I11.27P in the integral form by the Duhamel principle to obtain 



\\{x)-^P^'R{t)\\2 < |Kx)-'^e*^(^i)P,^ii?(0)||2 + /o ||(x)~'^e(*-«)-^(^i) 

X [P{z, z)R + PtS2{z, z) + P^^P^N2{R, z)]\\2ds. 

For the left hand side we claim that 

||(x)-^e*^(^^)p^P(0)||2 < c (1 + t)-i(||(x)'^P(0)||2 + |z(0)|2). 

Indeed, recall that 

R = R— ^ Rm,n 
m+n=2 



(11.28) 



(11.29) 
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with Rm,n defined in ( I9.15p . Therefore, displaying the time-dependent of R, A and z, 

m+n=2 

By ( I5.14P and the fact that the notation Rm,n is the summation of terms of order jzp we have 

||(x)-V^(^i)p,^ii?„,„(0)||2 <c|z(0)|2(l + t)-i. 
This together with the estimate 

||(x)-V^(^^)p^P(0)||2 <c(l + t)-i||(x)''i?(0)||2 

imphes f fTOO]) . 

Use ( I5.14P on the right hand side of Equation ( I11.28P to obtain 

< + t - srk\\{xyP{z, z)R\\2 + \\N2{R, z)||i + ||A^2(P, z)h)ds (11.30) 

+ /o II (x)-"e(*-^)^(^i)P,^i52(z, z)\\2ds. 

We estimate these terms in detail: 

(A) By the definition of S2iz,z) in Equation ( 18. 3p and Estimate ( I5.14p with d = 3 we have that 

I ||(a;)-'^e(*-^)^(^i)P,^^52(z,z)||2ds| < cj^{l + t - s)-i(To + syUsZ^ 



3 



< c{To + t)-^Z^. 

(B) By the definition of P{z, z) and the estimate of M2{z, z) in Equation ( 18.41) 

\\{xrP{z{slz{s))R{s)\\2 < C|Z||KX)--P(.)||2 

< c(ro + s)-iz7^3. 

Hence by (111. lip 

jl{l + t-s)-h{xYP{z{s),z{s))R\\2ds 

< cf^il + t- s)-l{To + s)~ldsT~^Zn3 

< c{To + t)-'^T~^>Zn3 

These together with ( I11.26P imply 

IKx)-^p,^^p||2 < c(i + t)-i[|Kxrp(o)||2 + kP(o)] 

+c(ro + ^)-iT(7^(Z7^l + Z7^3 + z^ + nl + nlUi) 

< c{To + t)-l[Tj\\{xrm\\2+Tj'\z\\0) 

+T^^ {zui + zus + z3 + 7^i + 7^i7^4)]. 

This implies ( I11.25P . 



□ 



Proposition 11.11. 

nl < \\R{o)fH2 + cTo^nl + z^TZi + z^nf + nlnl]. (11.31) 
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Before the proof we estimate the nonhnear terms. 
Lemma 11.12. 

\{{-A + i)p^jNiR,z),{-A + i)R)\ <c(^o + ^)-2[z27^l+7^27^|]. (11.32) 

Proof. As in Lemma 111.41 we decompose JN into the locahzed term Lac and the non-locahzed NonLoc 
defined in ( 18. 7p . The Locahzed part satisfies the estimate 

|((-A + l)Loc, (-A + 1)R)\ < c[|z|2 + Wixy'RWl]. 

By the definition of NonLoc in ( 18. 7p we obtain 

|((-A + l)NonLoc, (-A + 1)R)\ < c||(-A + l)i?||^||i?||^. 

This together with the definitions of estimating functions imphes ( lll.32p . □ 

Proof. By Equation ( I9.1ip , we have 

dt{{-A + l)R,{-A + l)R) 
= ((-A + l)iR, (-A + 1)R) + ((-A + 1)R, (-A + l)|i?) 



n=l 



with Kn, n = 1,2, 3, 4, defined as 

Ki := ((-A + 1)(L(A) + ^J)R, (-A + 1)R) + ((-A + 1)R, (-A + 1)(L(A) + ^J)R); 

K2 := A((-A + l)PcxR, (-A + 1)R) + A((-A + 1)R, (-A + l)PcxRy, 

Ks := -((-A + l)P^JN{R, z), (-A + l)R) - ((-A + l)R, (-A + l)P^JN{R, z)); 

■= ((-A + l)P^G, (-A + l)R) + ((-A + l)R, (-A + l)P^G). 

Recall the definition of the operator L[\) in ( 15. 4p . By the observation J* = — J and the fact that 
JL[\) is self-adjoint we cancel all the nonlocal terms in Ky. 

\Ki\ < c||(x)-^i?||^2 < c(To + t)-^nl 

By observing that |A| = 0(|zp) and PcxR is localized we have that 

\K2\<c\z{t)\^ \\{x)-'R{t)\\l,<c{To + t)-^Z\t)nl{t). 

By ( we have 

By the property of P^Q in ( 19. lip we have 

\Ki\ < c\z\'^\\{x)-''R\\h2 < c(To + ty^Z'^TZi. 
Collecting all the estimates above to obtain 

|^((-A + 1)R, (-A + l)R)\ < c{To + t)-^[nl + Z^TZi + Z^TZj + TZlnj]. 
After integrating the equation above from to i we have ( lll.3ip . □ 
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11.6 Estimate for Z{T) = max{To + t)2\z{t)\ 

Recall that by (FGR), eqn ^M), z*[Z{z,z) + Z*{z,z)]z > C\z\*. 
Proposition 11.13. There exists an order one constant m > 0, such that if 

m<To< \z{0)\-^, then (11.33) 

z{T) < 1 + z{T) ( z(r) + 7^?(^) + 7^|(^) + z(^)7^3(^) ) (11.34) 

Proof. By the equation ( I8.12P we have that 

^|z|2 = -z*[Z{z,z) + Z*{z,z)]z + Re{z Remainder {t)) (11.35) 
which can be transformed into a Riccati inequality 

dt\z{t)\^ < -C\z{t)\'^ + 2 \z{t)\ {Remainder {t)\. (11.36) 

By (ED, 

\z{t)\ \Remainder{t)\ < Z{T) ( Z {T) + nl{T) + nl{T) + Z {T)n^{T) ) , (11.37) 

where 5 = 2/5. We now use the following Lemma, proved in the Appendix 113.71 
Lemma 11.14. Suppose that z{t) is any function satisfying the equation 

dt\z{t)\^ < -\z{t)\^ + g{t), z{0) = zo (11.38) 

where g{t) is a function satisfying the estimate 

\9{t)\ < c# {To + ty-' (11.39) 

with constants c#, 6 > 0. Then, there exists K > 0, independent of Tq and c#, such that if c^T^^ is 
sufficiently small, then the function z{t) in ( \11. 3^) admits the bound 



\z{t)\ < ^ ' " ^^r>' (11.40) 

in+t)2 



1+ Kc# To 



where k = min{To, \zq\ 

We now apply this Lemma, choosing 

m < To < \z{0)\~^ (11-41) 
as in ( I11.33P where m is an order one positive constant. Then, 

Z{T) < 1 + ^ Z{T) {Z{T)+nl{T)+nl{T) + Z{T)n:i{T) ) 



□ 
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11.7 Closing the Estimates; Completing the Proof of the Main Theorem 17.11 

We seek to obtain T-independent bounds on TZj{T) and Z{T), defined in (jll.2p . This will be achieved 
by choosing the parameter Tq sufficiently large, appearing in the norm definitions, and the data -R(O) 
sufficiently small, with Tq and R{0) related in a manner to be specified. 
Define 

M{T) := 5^7^„(^), S := m{0)\\H2 + \\{xr R{0)h), (11.42) 

where, recall the definition of Tq = ( lll.ip . By the conditions on the datum ( 17. ip we have that 

7^4(0) is small and M(0) and Z{0) are bounded. 

Recah the estimates of 7^„, n = 1,2,3,4,5 and Z in Equations f [TlT3]1 . f [TlTH]) . r fTOHl) . f [TOT]) . 
(HOB) and ( flLMD . By plugging ( flLMD (flL3l|) . the estimate of Z and 7^4, into f nTTSll . ( flTTBD and 
( I11.21P we obtain 

M{T) < c{S + 1) + {R4{T) + T~^')P{M{T),Z{T)), 

Z{T) < 1 + T~^P{M{T),Z{T)), (11.43) 
7^1(^) < \\Rm\l,+T,'P{M{T),Z{T)) 

where P{x, y) > is a polynomial of variables x and y. Using an implicit-function-theorem type argument 
(see below) we have that if S and M(0) are bounded, then 

M{T) + Z{T) < i2{S) and 7^4 < 1, (11.44) 

where /i is a bounded function for S bounded. By the definitions of TZj{T) and Z{T) there exists some 
constant c such that 

\\{x)-''R{t)\\2, \\R{t)\\^<c{To + t)-\ \z{t)\<c{To + t)-^2 (11.45) 

which is Statement (B) in Theorem 17. li 

By the bound of Remainder in ( 18. 8p and the estimates ( I11.45P we have 

19 

\Remainder\ < c(To + t) 5" 

which together with ( I8.12P implies Statement (A) . 

The convergence of A comes from ( 18. 9p and the fact Remainder is integrable at 00. 

□ 

In the following we prove pi.43p implies ()11.44p by using implicit function theorem. For the other methods 
we refer to [171 [SH El HH [53l [SH |T7] . First we transform the inequalities by taking square root of the third 
equation of (I11.43P and plugging it into the first one, then there exists a fixed polynomial P{x, y) with 
positive coefficients such that 

M{T) < c{S + l) + mm\m+T;^')PiM{T),Z{T)), 

ZiT) < 1 + T~^P{M{T),Z{T)), (11.46) 
7^4(T) < ||i?(0)||H2+To"^P(M(T),Z(r)). 

In what follows we use this equation instead of (|11.43p . Define a vector function F^^siM, Z) as 

F,4M, Z) := {fP{M, Z), fP{M, Z), Ff}{M, Z)) 
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with 

F^^\m, Z) := c{S + l) + {6 + e)P{M, Z) 
Fp\M, Z):=l + eP{M, Z) and F^J -.= 5 + eP{M, Z). 
Immediately we can see that Mq = c(l + 5), Zq = 1 and i?o = is a solution to the equation 

iMo,Zo,Ro) = Fo,o{Mo,Zo,Ro). 

Define a closed set 

E := [0,2c(5 + l)] X [0,2] x [0,1]. 

Now we have the resuh. 

Lemma 11.15. There exists a 5o>0 such that if e,6 €z [0,6o] then 

(M, Z, R) = F,^s{M, Z, R) (11.47) 
has a unique solution in E, moreover for any continuous functions M, Z, TZ : — > M"*" satisfying 

(M(0),Z(0),7^(0)) < iM,Z,R) and {M (t) , Z (t) , R{t)) < F,^s{M{t), Z{t),n{t)) 
we have for any time t, 

{M{t), Z{t),n{t)) < (M, Z, R). (11.48) 
Proof. The proof of existence and uniqueness of the solution is not difficult by observing 

\\{dMF,4M, Z),dzF,4M, Z),dnF,4M, Z))\\ < c{6 + e) 

if (M, Z,TZ) G S. Hence by implicit function theorem we have that if c(e + S) < ^ there exists a unique 
solution to (|11.47|1 . 

We next prove pi.48p by contradiction. Suppose that pi.48|) fails at time t. Since (M(t), Z{t),n{t)) 
is continuous, there exists a time ti < t such that {M{ti), Z{ti),Tl{ti)) S E and (jll.48p does not hold. 
Without loss of generality we assume t = ti. Then by subtracting the inequality for {M{t), Z{t),Tl{t)) by 
(|11.47p we get 

M{t)-M < {S + e)[Ki{M{t)-M)+K2{Z{t)-Z)] 

Z{t) - Z, 7^(^) - ^ < eK^{M{t) -M) + eKi{Z{t) - Z) 

for some Kn, n = 1,2,3,4, depending on {M{t), Z{t),n{t)) and {M,Z,R). By the fact that 

{M,Z,R), {M{t),Z{t),n{t)) G E 

and P{x,y) is a polynomial with positive coefficient we have that Kn, n = 1,2,3,4 are positive and 
bounded. By these inequalities and the fact < e,(5 <C 1 we derive ()11.48p . This contradicts our 
assumption. Thus ()11.48p holds for any time t > 0. □ 
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12 Summary and Discussion 



We have extended the asymptotic stabihty / scattering theory of sohtary waves of the nonhnear Schrodinger 
/ Gross-Pitaevskii (NLS / GP) equation to the important case, where the hnearized dynamics about the 
Lyapunov stable has degenerate neutral modes. This is the prevalent case in situation, where the equation 
is invariant under a nontrivial symmetry. We construct a class of double-well potentials to which the theory 
applies. The current theory, as in all previous work on soliton scattering in systems with non-trivial neutral 
modes, requires a Fermi Golden Rule (FGR) non-degeneracy hypothesis. The analytical verification of 
this hypothesis for specific or generic NLS/GP systems is an open question. Numerical experiments for 
the time-dependent NLS/GP equations, in which decay rates of neutral modes is measured, are consistent 
with the typical validity (FGR) non-degeneracy hypothesis. 

We conclude by mentioning an interesting direction for further exploration. 
Semiclassical limits and higher order nonlinear Fermi Golden Rule: A problem of great interest is NLS / 
GP on M*^ in the semi-classical limit: 



The nonlinearity is taken to be focusing (attractive) but subcritical. Using the Lyapunov-Schmidt method 
it has been shown in [20^ [351 [1] that for h sufficiently small a soliton concentrated at a nondegenerate 
critical point of V can be constructed. The soliton, constructed in this manner, is soliton of the translation 
invariant nonlinear Schrodinger equation, scaled to be highly concentrated about the critical point of 
V. Therefore, the linearized operator JH^{X) is expected to have spectrum, quite closely related to the 
linearization about the translation invariant NLS soliton. If the soliton is concentrated near a minimum 
of V, then it is Lyapunov stable ^35j, and therefore the spectrum of JH^{X) is a subset of the imaginary 
axis. As we have seen for NLS/GP, there is a two-dimensional generalized eigenspace corresponding to an 
eigenvalue zero, h being small, implies that the 2 x d zero modes, associated with the translation symmetry 
(^(x, t) ^ 'ijj{x + xq, t)) and Galilean symmetry (^(x, t) giv-{x-vt) _ 2yf^ perturb to two complex 
conjugate eigenvalues, each degenerate, of multiplicity N = d. Although we expect semiclassical, highly 
localized solitons to be asymptotically stable and for the degenerate neutral modes to damp by resonant 
radiation damping, as elucidated in this article, we note that for h very small, the complex conjugate 
neutral modes of JH^{X) are very close to zero and the condition 2E{X) — A > 0, which is necessary 
(although not sufficient) for the Fermi Golden Rule resonance condition (FGR) to hold, fails. It remains 
an open question to derive the normal form, when resonance of discrete modes with the continuum occurs 
at some arbitrary order in the coupling parameter, g; recall /(|V'P)^ = see also the discussion in 

the introduction. For results in this direction, see \22\ I17j. 

13 Appendix 

13.1 A class of double-well potentials for which — A + y satisfies condition (Eigy) and 
L(A) satisfies (SA) and (ThreshA) 

In this section we find an example —A + y in a subspace of L^(]R^) satisfying condition (Eigy), motivated 
by the study of double well potentials. Define 



Observe that ^ is a self-closed subspace, i.e. if /i, f2 £ A then fi + /2, /1/2, A/i G A. Hence we can 
study (II. ip in the space L^(R^) and obtain all the results. The following is the main result 



idttp 



-Aijj + V{hx)7p- f{\ip\'^)ip, 0<h«l 
■tpoix) 



(12.1) 



A:={f :M.^ ^ C| f{-x) = f{x) for any x]. 
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Proposition 13.1. There exists a potential V such that the linear operator —A + V acting on the suhspace 
Ar\ L^(M^) has two eigenvalues cq < ei < with 2e\ > cq. cq is the lowest eigenvalue, the eigenvalue ei 
is degenerate with multiplicity 2. Moreover the operator 1 + {—A + iO)~^V : (x)^L^ — > (x)^L^ is invertible. 

If the nonlinearity f{x) = x, and if |A — |eo|| is sufficiently small and (j)^ is the ground state satisfying 

-A(/)^ + V(l)^ + X(p^ - {^^f = 

then we have the following results for the linearized operator L(X) defined in ()5.4p . 

Proposition 13.2. The operator L{\) satisfies the spectral conditions (SA) and (Threshx). 

The Proposition 113.11 is implied by Proposition 113.51 below. Proposition 113.21 will be proved at the 
end of this section. 

As proved in [2] (Theorem 1.1.4, Page 116) the operator —A — q6{x) has only one eigenfunction, i.e. 
the ground states, for any g > 0. By this observation we have 

Lemma 13.3. For any q> 0, there exists a constant A G (0, oo) such that the operators —A — qX ae ^ , 

, _ 3 aip 

— A — g^A 2e A both have only one eigenfunction in (pi, (1)2 ^ A respectively. 
To facilitate later discussions we define 

3 ^ 

W := gA~2 e ~ . 

We start with constructing a family of operators: Define Mi := (m,0, 0), M2 ■= (0, m, 0) and M3 := 
(0,0, m). And define Wm^ as 

Wm^x) ■■= ^iWix + Mk) + W{x - Mk)], k = l, 2,3. 

Now we have the following result. 

Lemma 13.4. If m is sufficiently large, then in the suhspace Af^Lp'i^) the operators —A — Wm,,, — A — 
^VFjv/fc; k = 1,2,3, each has only one eigenfunction. 

Proof. In that following we only prove the argument for —A — Wmi- The proof of the other cases are 
similar, hence omitted. 

First we have that if m is sufficiently large then 

((-A - WM,m- + Ml) + 0(- - Ml)], [,/.(• + Ml) + (/){■ - Ml)]) < 0. 

This implies that the operator —A — Wmi has at least one ground state. 

Secondly the min-max principle implies that any function / _L (/>(• + Mi), (j){- — Mi), 

((-A - WM,)f, f) = \ [((-A - W[- + M))f, f) + ((-A - W{- - M))f, /)] > 0. 

This together with the facts </>(• + Mi) - (/>(• - Mi) _L L2(]R3) n A and {</>(• - Mi), 0(- + Mi)} = {0(- - 
Ml) ± (/.(• + Ml)} yields that (((-A - WM,)f, /) > for any / G and / ± 0(- + Mi) + <j){- - Mi). 

Collecting what was proved we have that the operator —A — Wm\ has only one eigenfunction, i.e. the 
ground state. 

The proof is complete. □ 
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To prove the main result we have to define V. Define Vm as 

Vm:=^[WM,+WM,+WM,]. 

Now we have 

Proposition 13.5. There exists at least one m G [0, oo) such that —A — Vm has all the properties in 
Proposition \13.1\ 

Proof. To prepare for the proof we hst the following facts: 

(A) For any m G [0, oo) the operator —A — Vm has at most three eigenfunctions in ^ n L^. Recall 
in Lemma [HTi] we proved that if / -L </>(• + Mfc) + 0(- - M^.), A; = 1,2,3 and / G ^ n then 
((-A - WMk)f\ f) > 0. Consequently if / _L (/>(• + M^) + - M^), A; = 1, 2, 3 then 

((-A - Vm)f, f) = ^[((-A - WM,)f, f) + ((-A - WM,)f, f) + ((-A - WMs)f, /)] > 0. 

The min-max principle [37j implies that there are at most three eigenfunctions. 

(B) If m is sufficiently large, then in the space L'^ H A the operator —A — Vm has three eigenfunctions 
and two eigenvalues: one ground state and two (degenerate) neutral modes. The fact —A — Vm 
has three eigenfunctions comes from the min-max principle, the proof is similar to that of double 
well potential (see [291 132] ), hence is omitted here. The hard part is to prove the neutral mode is 
degenerate. Indeed, as m ^ oo the three eigenfunctions converges to some linear combination of 

+ Mk) + (pi- — Mk), k = 1,2,3 and the ground states converges to Ylt=i + ^k) + </>(• — Mk). 
Moreover the ground state must be simple and orthogonal to the neutral modes, i.e. the neutral 
mode can not be invariant under the permutation {xi, X2,x^) (a;„(i), x„(2)) 2;„(3)) while the operator 
—A — Vm is. This enable us to obtain one of the neutral modes by permuting any one of two, hence 
the eigenvalue of neutral modes must be the same. 

(C) When m = 0, —A — Vm has only one eigenfunction, the ground state. This is obvious by the fact 
Vm = W when m = 0. 

(D) For any m > 0, —A — Vm has at least one eigenfunction with eigenvalue less than some — cq < 0. 
Let 02 be the normalized ground state of —A — \Wm2 with eigenvalue — cq < 0. Then we have 
((—A — Vm)4>2-, 4'2) < —Co by the facts 02 > and W > 0. By the min-max principle —A — Vm has 
one ground state. 

The definition of W implies that (—A + k)^^W{- + z) is analytic in z if G C\M+. By ^37j we have 
that the eigenvalues are analytic functions of z in a suitable subset of C. Since the eigenvalue of the 
neutral modes is degenerate for sufficiently large m (see (B)), it is degenerate for any m before the neutral 
modes disappear into the essential spectrum. Hence there exists at least one m such that —A — Vm has one 
eigenvalue less than — cq (defined in (D)) and two degenerate neutral modes with one eigenvalue sufficiently 
close to the essential spectrum (see (A), (C)). 

In the final step we find m and q such that the operator 1 + (—A + iO)~^Vm '■ (x)^L^ 
is invertible. Recall that Vm = qV2{m) with V2(m) independent of q by its definition. For a fixed qo 
we proved that there exists at least one m = mo such that the eigenvalues of —A — qoV2{mQ) have the 
desired properties. Now we consider a family of operators X{g) := 1 — q{—A + iO)~^l^(mo) which is 
analytic in q and the operator (?(— A + iO)~^V2{mQ) : {x)'^L? — > (x)^L^ is compact. By [37| the operators 
X{q) : (x)^L^ — > (x)^L^ are either invertible everywhere (i.e. no threshold resonance) except discrete 
points, or are not invertible anywhere. The first case holds because the operator is invertible when q = 0. 
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Now we consider —A — qV2{mQ) with q £ [qo — e, qo + e]. Choose e sufficiently small such that for every 
q the operator —A — gV2(?TT-o) has at least three eigenvectors. On the other hand by what we proved above 
it has at most three eigenvectors and the second eigenvalue must be degenerate. Now by applying the fact 
that I- q{-A + iO)-^V2 (mo) is not invertible only at discrete points we obtain the desired result. 

The proof is complete. □ 

Now we prove Proposition 113.21 
Proof of Proposition ll3.2I The fact L{X) has no resonances at itiA is due to the facts l + {—A + iO)~^V 
is invertible and |A — |eo|| is small, see Proposition 15.31 

In the next we prove the neutral mode is degenerate. Recall that the potential we constructed is of the 
form V = for some mo- For each m > there are A = A^ and (j)^ = cj)^''^ satisfying the equation 

-AcP^'"" + A„(/.^''" + - i^^n^ = 

with Am and (pm analytic in m in some proper neighborhood of positive real axis. 

Recall that when m is sufficiently large the neutral modes of —A + Vm can be generated by permuting 
one of them, hence the neutral modes of L{\) = L(X, m) are degenerate when ra is large. Moreover the 
eigenvalues of L(A, m) are analytic in m, thus the neutral modes must be degenerate. 

The proof is complete. 

□ 

13.2 Fermi Golden Rule - Proof of Theorem W7\\ 

The proof of Theorem 16.11 uses the following 

Proposition 13.6. Given smooth functions T ^ Q : M"' — > C^, there exists T = {■^i,.^2) o.nd Q = (^1,^2) 
(see definitions below), such that 

- Re ((L(A) + 2iE{X) - oy^P, iJPc G) = tt 6{-A - {2E{X) - X) )J^2, Q2 ) (13.1) 

The proposition will be proved later. Proof of Theorem I6.lt We use Proposition 113.61 with T = Gk 
and G = Gi, ^ = Gk and Q = Gi. By (jl3.ip we have 

Tk,i = vr ( 5{-A - {2E{X) - X) )Gi,2, Gk,2 ) ■ (13.2) 

To see that Tj^^i is non-negative, observe that for any s G C''^, we have 

N 

s*Ts = J2 ^k,isksi = TT ( 5(-A-(2^(A)-A) )g,g) > 0, (13.3) 

k,l=l 

where G = J2k=i Gk,2- 

For the second statement we only sketch the proof. Recall the transformation of L{X) in (15. 7p . Then 
for any 2x1 vector functions F and G we have 

{{L{X) + 2iE{X) -Oy^PcF,PcG) = -i{{H + 2E{X) +iO)-^A*PeF,A*PcG) 

= -i{K{X){Ho + Vi + 2E{X) + iQ)-^A*PcF, A*PcG) 
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with K{X) is the operator defined as [1 + Ksmaii]~^ with Kgmaii ■= {Hq + Vi + 2E{X) + iO)~^Vsmaii- The 
operator {Hq + Vi + 2E{X) + iO)~^ is well defined by the fact —A — 2E{X) eo — 2(ei — eo) is not an 
eigenvalue of —A + V and the operator —A + V has no embedded eigenvalues in the essential spectrum. 
Since the operator K small '■ {x)'^L°^ — > {x)'^L°° has a small norm and continuous in A we have 

oo 

[1 + K small] ^ = y^^j — Ksmall)"' 
n=0 

is continuous in A. This together with the fact {Hq + Vi + 2E{X) + iO)~^ A* PcF G {x)'^L°° is continuous in 
A implies that {{L{X) + 2iE{X) - 0)-^PcF, PcG) is continuous in A. 
The proof is complete. 

□ 

Proof of Proposition 113. 6t The entries of T are expressions of the form 

- Re ((L(A) + 2iE{X) - O^^Pc iJPc Q) , (13.4) 
which we now proceed to simplify. Recall L{X) is of the form 

where Vi and V2 are real- valued and exponentially decaying as |x| tends to infinity. Introduce the unitary 
matrix 

J_ / 1 i 



Note that 

L(A) = iU as n{X) U* , H* = H (13.6) 



where 



n := no+ V 

Ho ■■= (-A + A) / (13.7) 



V 



0-3 := 



Vi - V2 -i{Vi + V2) 
i{Vi + V2) V1-V2 

-1 ; 

We now use the unitary transformation, U, to obtain an expression in terms of the operator (J3H: 

- ((L(A) + 2iEiX) - oy^Pc iJ Pc Q) 
= - (^{iU [asHiX) + 2E{X) + iO] U* )"^ Pc T, iJ P^ 

= (asnix) + 2Eix) + ioy^ u*Pc T, u*j Pc q) 

= {(JsHiX) + 2E{X) + iOy^ U*Pc {U*J U) U*Pc 

iasTiiX) + 2E{X) + iOy^ U* Pc (J3 U*Pc q) , (13.^ 
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where we have used that U*JU = ia^. 

Next, we introduce Pcicr^TC), the projection onto the continuous spectral part of a^H and wave operators 
: L2 ^ Pc{a3n)) L2 and Z : Pc{a3H)L'^ L^, see HI], which satisfy 

Pci^sny <y^ = Pc{cT^n), 

W* as = 0-3 Z, Z* as = asW 

Z asT-i = asTio Z (13.9) 

Now we use the wave operators W and Z to transform the previous expression into one, in terms of the 
"free operator" (T3(— A + A). First, note that U*Pc T hes in the range of Pda-sTC), and therefore there exists 
jr = {Ti,T2f such that = U*Pc J=. Similarly, there exists Q = (^1,^2)^, such that WQ = U*Pc Q. 
Substitution into the final expression in (jl3.8j) and use of the properties (jl3.9p we have: 

i ( (cT3W(A) + 2E{\) + zO)-^ U*Pc J", as U*Pc q) 
= i ( iasniX) + 2E{X)+iO)-^ WT,as WQ^ 
= i ^ iasniX) + 2E{X)+iO)-^ W^,Z* as q) 
= i i^Z {asUiX) + 2E{X) + iO)"^ WT, as g) 
= i (^{as{-A + X) + 2E{X)+i0y^ ZWT, as q) 

= i ^ {as{-A + X) + 2E{X)+i0y^ T, as q) (13.10) 
Referring back to (113. 4p . we recall that we are interested in the real part of this expression. 

-Rei {as{-A + A) + 2E{X) + iOy^ as q) 

= Im l^{as{-A + X) + 2E{X) + iQy^ T, as O) 

/ / (_A + A + 2i?(A)+iO)-i \ \ 

= ^"^ U -(-A + A - 2E{X) - ^0y^) ) / 

- Im/ ( (-A + A + 2i?(A)+zO)-i W c\ 

- ^"^[[o {-A + x-2E{x)-ioy') ) ^ / 

= Im (-A - (2^(A) - A) - i0)~ij^2, ^2 ) 

= TT 6i-A - (2^(A) - A) )^2,G2 ) . (13.11) 

The last equality uses that < 2£'(A) — A G (Tc(— A), — 2ii^(A) — A ^ cTc(— A) and the distributional (Plemelj) 
identity: 

Im {x — iO)^^ = lim Im (x — ie)~^ = vr S{x) (13.12) 

Summarizing, we have shown 

- Re ((L(A) + 2i^(A) - Oy^P^ iJPc G) = n 6{-A - {2E{X) - A) )J^2, ^2 ) • (13.13) 
This completes the proof of Proposition 113.61 
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13.3 FGR for symmetric potentials 

In this section we derive the simpler form of the FGR matrix and condition for positivity in the case where 
the potential V{x) is a function of \x\. In fact, it is proved in Lemma 15.21 that if the potential V, hence 
(p^, is spherically symmetric, then the functions r/„ satisfy 

= Vn = ^r?(|x|) (13.14) 

\x\ \x\ 

for some functions and r?(|x|). By the assumptions on V, cj)^^ ^k, Vk, k = 1,2, ■ ■ ■, N = d we have 

Gkiz,x) = Xkiz ■ x)G{\x\) 

for some radial vector function G(|a;|). 

Before stating the results we define two constants 

Re Z^^'^'^ = -Re {iL{X) + 2iE{X) -Oy^Pc xj G{\x\),i J xf G{\x\)) 

Re Zp'^'^ = -Re {{L{\) + 2iE{X)-{))-^PcXiX2G{\x\),i J xiX2G{\x\)) (13.15) 

Proposition 13.7. (i) Suppose that V , ^n, ??n satisfy the conditions above. Then the assumption (FGR) 
holds provided that 

ReZ^^'^^ > 0, ReZl^^'^^ > 0. (13.16) 
(a) From. Proposition \13.6\ it follows that 

i?eZi^'^^ > 0, i?eZi^'^^ > 0. (13.17) 

and, generically, that strict positively il3.16\) holds. 

Proof: For any vectors s,f3,z S C^, we define 

Q{s, p; z) ■= -Re ((L(A) + 2iE{\) - Q)~^Pc{z ■ x){s ■ x)G{\x\),iJ{z ■ x){(3 ■ x)G{\x\)) , (13.18) 

Note that ^ 

Q{s, s; z) = -s* [Z{z,z) + Z*{z,z)\s = Re s*Z{z,z)s 

Therefore, verifying (FGR) is equivalent to checking that there is a constant C > 0, for which 

Qis,s;z) > C |sp s,z£&. (13.19) 

To simplify Q{s, s; z), first note that since operator L{X) and G(|x|) are invariant under transformations 
x T*x, where T is unitary, the value of Q{s, (3; z) is unchanged when replacing x by T*x. 
Therefore, 

Q{s,P;z) = Q{Ts,TP;Tz) (13.20) 
Now choose T to be a unitary matrix, such that 

Tz = \z\ ei = \z\ (1,0, ...,0)'^. 

With this choice of T, we have by (jl3.20p with (3 = s, 

Q{s,s;z) = -Re{{L{\) + 2iE{\) - {))-'^P^\z\xi{Ts ■ x)G{\x\),i\z\xi{Ts ■ x)JG{\x\)) . (13.21) 
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The following argument will show that Q{s,s;z) > C \Ts\ \z\^ = C \s\^ jzp, the latter holding since 
T is unitary. Therefore, without any loss of generality consider ()13.2ip with T set equal to the identity. 
Explicitly writing out the inner products and using bilinearity and symmetry, we have 



Re ^ {{L{\) + 2iE{\)-{))-^PcXiXpG{\x\),ixiXyJG{\x\)) SpTg 



d 

|2 



\Z\ 

p,q=l 
d 

Re ((L(A) + 2iE{X) - Oy^Pc xiXp G{\x\),i xiXp JG{\x\ 
p=i 

\z\^ Re {{L{X) + 2iE{X) - or'^Pc xl G{\x\),i xj JG{\x\)) \sif 



z 



Re ((L(A) + 2iE{X) - oy^Pc xiX2 G{\x\),i xiX2 JG{\x\)) ^ | 

q=2 

q=2 



> mm{Re Z^^'^\ Re Z^^'^^ } |zp = C \z\'^ > 0. 
The proof is complete. 

13.4 Choice of Basis for the Degenerate Subspace: Proof of Proposition 15.41 

In the proof of the proposition we need the following lemma. 

Lemma 13.8. If u = { ' / is an eigenf unction of L(X) with eigenvalue iE{X), E{X) > 0, then 

\ J 

{ui,U2)>0. (13.22) 

Proof. The fact L(X)u = iE{X)u yields 

L_{X)u2 = E{X)ui, L+{X)ui = E{X)u2. (13.23) 

Therefore ^ 

(ni,n2) = -g^(L_(A)n2,n2). 

Equation ( 113.220 follows from the two claims that L-{X) is a positive definite self-adjoint operator on 
the space {v\v _L cp^} and U2 span{(j)^}. The first fact is well known (see e.g. [57]). We prove the second 
by contradiction. Suppose that U2 = c(f>^ for some constant c, then we have L-{X)u2 = 0, which together 
with ( 113. 23p and the fact E{X) ^ implies ui = U2 = 0, i.e. n = 0. This contradicts to the fact u ^ 0. 
Thus U2 span{(j)^}. □ 



Proof of Proposition l5.4I We start the proof by constructing N independent vectors n„ G span{vi, V2, ■ 

1 

i 



•, f at}, n = 1, 2, • • •, such that the vector ( ^ ^ ] is real 
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Suppose that = I (n) 1 • Then the definition of L(A) in ( 15. 4p implies I ^ («) 1 I . ^ (n) 



ilmv2 I \ —iRev\ 



are also eigenfunctions of L{\) with eigenvalues iE{\). This together with the fact 

,n = l,2,---,A^} = K, n = l,2--,A^} 



enables us to choose independent eigenfunctions for iE{\): Un, n = 1,2, ■ ■ N, such that ^ ^ ^ j Un 
are real vectors. 

Using ( 113. 22p and a standard Gram-Schmidt procedure in linear algebra, one can find N pairs of real 
functions rj^), n = 1,2 • ••,A^, such that span{( ^" | , n = 1,2 • ■■,N} = span{vn,n = 1,2,--,A^} 

and {in-,rira) = 5n,m- 

We now turn to the verification of ( I5.10p . The observations 
L_(A)-L+(A) = 2f'U^f]{ct>^f, 

L_(A)r/„ = E[\)in, andL+(A)e„ = ^(A)?7„, n = 1, 2, • • •, A^, 

yield 

! fU^?]{4>''?{Urin-inrim)dx 
= 0. 

Finally, ( 15. lip is seen as follows: 

{^\^n) = -^(0\L_(A)r?„) = -J^(L_(A)<A\r?„) = 0, 

13.5 The identity P^{JHyj = JPc{JH) 

Proposition 13.9. L = JH,H = H* =^ [PciL)]*J = JPc{L) . 

Proof: Represent Pc{L) as a Riesz projection 

PJL) = — lizi - JPy^dz, (13.24) 
2m J 

where the integration is counter-clockwise, moreover the essential spectrum of L is [iA,ioo) U {—ico, 
The spectrum associated with the upper branch, [i\,ioo), is given by 

P+{JH) = —[A-B], where (13.25) 
27r 

/•oo 

A = J (iT + O-JH)-^ dT 

roo 

j {iT -0- JH)-^ dr. 



B 

We claim that 



A*J = -JB,B*J = -JA. (13.26) 
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This implies [P+{JH)]*J = JP+{JH). Similarly, [Pc{JH)]*J = JPc{JH). 

To complete the proof of Proposition 113.91 we now prove (|13.26p . By direct computation, using 
J* = — J, we have 



j-OO 

^* = J i-ir + + HJ)-^ dT 



Therefore, 



A*J = / {JirJ - JOJ - JJHJ)-^ dT J = {-J){iT -0- JH)-^{-J) dT J = -JB. (13.27) 
Jx Jx 

thus proving the first identity in (jl3.26p . B*J = — J A can be proved similarly. 
13.6 Time Convolution Lemmas: Proof of Proposition 111.21 

Proof: In what fohows we only prove the case o" = 1 of ( lll.lip . the other cases and ( I11.12p are similar. 

w X /■* 1 1 , 1 /"^ 1 , If*! 
I(t) := / ^ ds < ^ / ds + r / 

Jo {l+t-s)-2 n + s (1 + 1)170 To + s To + I (i+t_s) 

log(l + 27b 



^ + 



On the other hand, we also have 



Thus, 







(l + |)i ^0 + 1 

* 1 1,2 
3" ds < — 

{l + t-s)2 To + s To 

1 1 



I(t) < ci min 



To ' 1 + t 



We now claim that for some constant c > 0, that I{t) < c (To + t) ^. It sufficies to find a constant c, 
independent of Tq and t, such that 

1 1 

m{t) := (To + 1) min{ — , — — } < c 

io i + r 

If t is such that the above min is Tq^, then ro""*^ < (l + t)'^^ or t < To — 1. Then, m{t) < (2To — 1)^0^^ < |. 
If t is such that the above min is (1 + then t > To — 1. Therefore, m{t) < (2To — ^)Tq^ since m{t) is 
decreasing with t. Since T > 2, m{t) < 3/2. This completes the proof. 

13.7 Bounds on Solutions to a Weakly Perturbed ODE: Proof of Lemma 111.141 

Proof. Let /3 denote the solution to the differential equation 

\Pp\\0) = |z(0)|2 - p, p>0. 

Since 

dt{\zit)\'-\p,{t)\') < + 
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with the initial condition 

|z(0)|2-|/3,(0)|2 = p>Q. 
Thus, < |/?p(t)p for all t > and letting p tend to zero, we have 

so it suffices to prove the bound: 

|/3(t)| < (l + iC c# Tq-^) + K = min{ro,|u;or'} (13.28) 
where f3{t) solves the initial value problem 

\m? = kop. (13.29) 

The proof of (jll.40p for (3 is divided into two cases: \wq\ > Tq and |?x;o| < 

First consider case (i): |w;o| > Tq . By local existence for the initial value problem (jl3.29p . we have 
that for some ti > 0, 

< m% tG[0,ti]. (13.30) 

Then, using the assumed bound on g{t), (|11.39p . we have 

c# _ c# 1 



\9{t)\ < 



(ro + t)2+^ (To + t)2 {n + tf 



where ci^ := 2^ c^. It follows from (jl3.29p that 

5d/3(i)P < -(l-ci#V)|/3(t)|4 

or 

-2 \ / 1 „ rp — (5 



dt\mr > { i-ci#ro 

Integration over the interval [0,t], t < ti, yields 



w{t)\ < ; ' (13.31) 



1 + C2# 



where C2# ~ ci# ~ c#, and we use that c# Tg is sufficiently small. Now set k = mm{\wo\ ^,To} and we 
have ^ 

\z{t)\ < m)\ < ^+''^#'^0 o<t<t. 

Now suppose that [0, H) denotes the maximal subset of IR+, on which the upper bound in (jl3.3ip holds. If 

_ 1 

H < cxD, then by continuity and the assumption that \wq\ >Tq ^, we have 

m.^M _ 1 + C2# V > 1 > 3 1 



(|y;o|~2+H)2 (|w;o|-2+H)2 4 (ro + H)2 

(13.32) 



52 



implying (see (|13.30p that the above argument can be apphed beyond t = H, contradicting its maximaUty. 

— 1/2 

Finally, we consider case (ii): |/?o| < Tq 

— 1/2 

Denote by Pi{t) the solution to ( lll.38|) with the initial condition /?i(0) =Tq .As shown in the previous 
case \Pi{t)\ <{1 + K c#T^^){Tq + t)^^/^^ Observing that 

dt{W-W) = -(|/3p + |/3i|2)(|/?|2-|/3i|2), 

|/3op-|/3i(0)P < 

we have that for any time t, < This together with the estimate of /3i completes the proof 

of the second case. □ 
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